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ABSTRACT 

Carmichael  C12  has  proposed  a  direct  injection, 
temperature  ignition  <TI),  reci procati ng  engine  design  which 
divides  the  conventional  -four  stroke  cycle  functions  between 
a  compression  cylinder  and  an  expansion  cylinder, 
interconnected  by  a  static  regenerative  heat  exchanger.  The 
prediction  of  heat  transfer  loss  in  the  expansion  cylinder  is 
required  to  predict  the  performance  of  the  new  cycle. 

Woschni's  correlations  C93,  developed  for  direct 
injection,  quiescent,  four  stroke,  compression  ignition  (Cl) 
engines,  are  used  to  predict  uncorrected  expansion  cylinder 
heat  transfer  loss.  Laws  of  similarity  allow,  in  principle, 
extrapolation  outside  the  experimental  data  range  for  forced 
convection  heat  transfer  only.  General  applicability  of 
Woschni's  correlations  is  established  by  examining  the 
expansion  cylinder  quiescent  condition,  using  predicted  mean 
gas  velocity  magnitudes  and  fully  developed,  turbulent  flow, 
forced  convection  models.  Predicted  radiant  peak  flame 
emissivity  and  temperature  are  found  and  compared  with 
nominal  peak  values  observed  in  conventional  Cl  engine  cycles 
to  establish  a  minimum  radiation  flux  reduction  ratio.  A 
mean  reduction  ratio  for  forced  convection  flux  is  also 
found.  A  corrected,  conservative  prediction  of  heat  transfer 
loss  in  the  expansion  cylinder  is  then  established. 
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NOMENCLATURE 


heat  transfer  surface  area  Cft^D 

cylinder  heat  transfer  surface  area  [ftA] 

regenerator  port  exit  cross-secti onal  area 
Cf  t^3 

crank  radius  length  Cftll 
bore  diameter  CftD 

fuel  consumption  rate  tlbm/hr] 

Boltzmann  number;  defined  by  equation  (61) 

friction  coefficient;  defined  by  equation  (11) 

specific  heat  at  constant  pressure 
CBTU/lbm-0R3 

specific  heat  at  constant  volume 
CBTU/lbm-°R3 

Woschni  correlation  constants;  defined  by 
equations  18a  -  13e 

speed  of  light  =  2.997902  x  1010  cm/s ec 

Boltzmann's  constant  =  1.38026  x  10-lto  erg/°K 

Planck's  constant  =  6.623777  x  10-ji7  erg-sec 

soot  formation  rate  coefficient 
C 1 bm/sec-f t°-atmD 

diameter  Cft3 

radiant  emissive  power  CBTU/hr — ft^  or 
BTU/sec-f t^] 

radiant  monochromatic  emissive  power 
Cergs/sec-cm^-?*  (cm)  3 

soot  combustion  activation  energy 
CBTU/mol e <lbm) 3 

radiant  geometric  view  factor 

heat  transfer  coefficient  CBTU/hr-f ti-°R] 


heat  transfer  coefficient  across  pipe^ 
diameter  D;  turbulent  flow  CBTU/hr-f t"*— °R3 


nrpm 


total  cylinder  heat  transfer  coefficient; 
turbulent  flow  CBTU/hr-f  ^-“Rl 

heat  transfer  coefficient  across  length  1  of 
flat  plate;  turbulent  flow  CBTU/hr-f t2-°R3 

local  heat  transfer  coefficient  at  distance 
from  leading  edqe  of  flat  plate;  turbulent 
flow  CBTU/hr-f t^-°R3 

Colburn  factor;  defined  by  equation  (9) 

thermal  conductivity  CBTU/hr-f t-° R3 

apparent  flame  absorption  coefficient  per 
unit  length  Cft'1] 

apparent  gas  absorption  coefficient  per  unit 
length  Cft-13 

apparent  flame  absorption  coefficient  per 
unit  length-unit  pressure  Cft-atm3-1 

monochromatic  absorption  coefficient  per  unit 
length  Cft-MJ  1 

piston  stoke  Cft 3 

cylinder  gas  volume  height  Cft3 

flame  luminosity 

flat  plate  length  Cft3 

Re  characteristic  length  Cft3 

effective  radiant  mean  beam  length  Cft3 

connecting  rod  length  Cft3 

charge  mass  Clbm3 

charge  mass  at  ignition  start  Clbm3 
charge  mass  flow  rate  Clbm/sec3 
cylinder  cycle  frequency  CRPM3 
Nusselt  number;  defined  by  equation  (3) 
cylinder  pressure  Catm3 


a 


cylinder  pressure  at  ignition  start  Eatm3 


PQ  cylinder  motoring  pressure  Eatm3 

Pr  Prandtl  number;  de-fined  by  equation  (5) 

Pu  unburnt  fuel  partial  pressure  in  the  scot 

formation  zone  Eatm3 

Q  heat  transfer  EBTU3 

Q  heat  transfer  rate  EBTU/sec3 

q  heat  transfer  flux  CBTU/sec-f t^l 

q^  monochromatic  radiant  heat  transfer  flux 

CBTU/sec-f  t"--,\3 

Rc  crank  ratio  =  lr/a 

RQ  universal  gas  constant  EBTU/mol e ( 1 bm) R3 

R  specific  gas  constant  =  53.90  lbf-f t/lbm-°R 

Re  Reynolds  number;  defined  by  equation  (4) 

ReD  Reynolds  number  for  pipe  flow;  lc  =  D 

R®x  local  Reynolds  number  for  flat  plate  at 

distance  x;  lc  =  x 

rc  cylinder  compression  ratio 

Sp  piston  speed  Eft/sec3 

A  soot  concentration  Elbrn/ft^l 

SFC  specific  fuel  consumption  Elbm  fuel/HF'-hr3 

St  Stanton  number;  defined  by  equation  (10) 

s  straight  line  distance  between  cylinder  crank 

shaft  and  piston  pin  Eft] 

T  temperature  E°R3 

Tg  cylinder  bulk  gas  temperature  E°R3 

Tgc  cylinder  bulk  gas  temperature  at  ignition 

start  C°RJ 

TQ  cylinder  motoring  bulk  gas  temperature  E°R3 
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local  temperature  of  the  soot  -formation  zone 
C  0  R  3 


u 

Tw 

t 

V 

Vc 

Vc 

Vu 

U 

Ui 

Uexp 

Uexh 


ymbol  s 
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a 
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heat  transfer  surface  wall  temperature  C°R3 

unit  time  Csec3 

cylinder  gas  volume  Eft^3 

cylinder  clearance  volume  at  piston  TDC  Eff'3 

cylinder  volume  at  ignition  start  (ft"') 

volume  of  the  soot  formation  zone  Eft'~'3 

bulk  fluid  (gas)  velocity  Eft/sec3 

impulse  charge  velocity  exiting  regenerator 
port  Eft/sec3 

bulk  gas  velocity  at  the  cylinder  walls  in 
the  expansion  stroke  Eft/s ec3 

bulk  gas  velocity  at  the  cylinder  walls  in 
the  exhaust  stroke  Eft/sec3 

angular  gas  velocity  at  radius  r  Cft/secD 

radiant  path  length  Cft3 


charge  equivalence  ratio 

unburnt  fuel  equivalence  ratio  in  the  soot 
formation  zone 

crank  angle  C°  degrees] 

forced  flow  boundary  layer  Eft 3 

electromagnetic  wavelength  Ccm  or  >un3 

dynamic  viscosity  Elbm/ft-sec3 

density  Clbrn/ft0’] 

kinematic  viscosity  (momentum  diffusivity)  = 
}i/p  Cft^/secl 

thermal  diffusivity  =  K/CpP  Eft^/sec3 
ratio  of  specific  heats  =  Cp/Cv 

lO 


A  heat 


apparent  shear  stress  acting  on  the  fluid  at 
the  wall  Clbf/ft2] 

emissivity  coefficient;  grey  body 

apparent  flame  emissivity  coefficient;  grey 
body 

flame  emissivity  coefficient;  grey  body 

gas  emissivity  coefficient;  grey  body 

emissivity  coefficient  of  infinitely  thick 
radiating  flame;  grey  body 

wall  emissivity  coefficient;  grey  body 

transmittance 

absorptance 

monochromatic  absorptance 
ref  1 ectance 
wall  reflectance 

—9 

Stef an-Bo^tzmann  constant  =  1.713  >•;  10 
BTU/hr— f t2— °R4 

incremental  change  or  difference 

specific  heat  release  per  cycle 
CBTU/lbm-chargel 

heat  transfer  loss/heat  release  per  cycle 

nominal  Qp./Qj  per  cycle  in  conventional 
four  stroke,  high  rc,  direct  injection  Cl 
engines 
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Subscr ipts 

cp 

combustion  products 

FC 

■forced  convection  heat  tran( 

NCE 

New  Cycle  Engine 

R 

radiation  heat  transfer 

T 

total  heat  transfer 

W 

Woschni  test  engine 

£ 

crank  angle  degrees 

Superscripts 

- 

mean  value 

* 

corrected  value 

Abbrevi ations 

BDC 

battom-dead-center 

TDC 

top-dead -center 

BTDC 

before  TDC 

C/H 

carbon/hydrogen  mass  ratio 

Cl 

compression  ignition 

SI 

spark  ignition 

T I 

temperature  ignition 
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CHAPTER  ONE 


INTRODUCTION 

1 .  SCOPE 

Conventional  compression  ignition  (Cl)  reci procati  ng 
engines  have  -found  significant  application  in  commercial, 
industrial,  and  in  land  and  marine  tr anspor tati on  systems. 
Recent  emphasis  on  fuel  costs  and  aval labi 1 i ty,  in 
conjunction  with  improved  materials  engineering,  has  prompted 
consi derati on  of  design  alternatives  which  offer  increased 
thermal  efficiency  and  reduction  in  operating  costs.  A  New 
Cycle  Engine  design  has  been  proposed  by  Carmichael  Cl  3  in 
which  the  functions  of  a  conventional  four  stroke 

reciprocating  cylinder  are  divided  between  a  compression 

% 

cylinder  for  charge  induction  and  compression,  an  expansion 
cylinder  for  charge  combustion-expansion  and  exhaust, 
interconnected  by  a  static  regenerative  heat  exchanger. 

The  attracti veness  of  such  an  engine  is  possible 
because  improved  fabrication  techniques  presently  allow 
construction  of  strong,  thermal  resistant,  extremely  thin- 
walled  ceramic  matrices  possessing  high  heat  transfer 
properties.  A  cyclic  interchange  of  compressed  charge  into 
and  expired  exhaust  gas  out  of  the  expansion  cylinder  via  the 
regenerator  establishes  an  effective  medium  for  heat  transfer 
and  allows  use  of  exhaust  heat  otherwise  lost  to  the  system. 
The  temperature  of  charge  air  heated  by  adiabatic  compression 
is  boasted  from  the  regenerator  heat  exchange,  a  direct 
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consequence  of  which  is  significant  reduction  in  cylinder 
compression  ratio,  r  ,  required  to  achieve  optimum 
eff iciency.  A  computer-generated  thermodynamic  model  of  this 
engine  in  preliminary  design  indicates  that  improved  thermal 
efficiency  may  be  achieved.  An  important  element  in  the 
cycle  analysis  is  accurate  estimation  of  heat  transferred 
from  the  gas  medium  in  the  expansion  cylinder  to  its 
enclosing  surface  area,  and  ultimately  to  the  coolant.  The 
magnitude  of  this  heat  transfer  reduces  both  heat 
availability  to  the  regenerator  and  the  thermal  efficiency. 

The  purpose  of  this  thesis  is  as  follows: 

1.  To  evaluate  the  applicability  of  existing 
heat  transfer  correlations  developed  to 
predict  total  mean  instantaneous  heat 
transfer  coefficient,  he,  in  four  stroke 
direct  injection,  quiescent,  Cl  engines 
for  prediction  of  total  heat  transfer  in 
the  New  Cycle  Engine  expansion  cylinder. 

2.  To  evaluate  and  compare  the  magnitude  of 
the  reduction  ratio  in  both  the  forced 
convective  and  the  radiant  heat  transfer 
from  the  high  peak  combustion  pressure 
conditions  inherent  in  high  rc,  Cl 
engines  to  low  peak  combustion  pressure 
conditions  inherent  in  the  low  rc,  New 
Cycle  Engine  design. 
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To  use  the  reduction  ratios  to  correct 


existing  heat  transfer  correl ati ons,  1 f 
applicable,  and  establish  a  correction 
allowing  better  prediction  o-f  heat 
transfer  loss  in  the  New  Cycle  Engine 
expansion  cylinder. 

2.  NEW  CYCLE  ENGINE  DESIGN  BACKGROUND 

A  vertical  cross-secti on  o-f  the  New  Cycle  Engine  design 
is  provided  in  Figure  1-1.  For  simplicity,  each  cylinder  is 
considered  o-f  simple  geometry,  with  straight  -forward  walls 
and  -flat  cylinder  head/piston  crown  configuration;  bore 
diameter,  B,  equals  stroke  length,  L.  Gas  velocities  are 
boosted  by  piston  speed,  induction  suction,  charge  impulse 
exiting  the  regenerator  port,  and  combustion  turbulence  only; 
quiescent  flow  should  predominate  with  no  deliberate  design 
feature  incorporated  to  accelerate  any  rotational  velocity 
incurred  by  off-set  regenerator  port  impulse.  The  cylinder 
crank  shafts  will  be  mechanically  linked  to  enable  the 
compression  piston  to  cycle  at  twice  the  speed  of  the 
expansion  piston,  alternating  between  an  induction- 
compression  and  a  free  motoring  stroke.  Proper  crank  angle 
phase  shift  will  exist  to  sequence  the  compression  piston 
motion  at  the  start  of  expansion  piston  down-stroke  to 
optimize  charge  transfer. 

Correct  value  timing  is  imperative  and  is  sequenced 
through  one  charge  cycle  as  follows; 
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1.  Compression  cylinder  induction  s' roke: 
induction  inlet  valve  1  opens; 
regenerator  transfer  valve  2  is 
c 1 osed . 

2.  Compression  cylinder  compression 
stroke;  valve  1  closes;  valve  2 
remains  closed. 

3.  Expansion  cylinder  combustion- 

expansion  stroke:  regenerator  exhaust 
value  3  closes;  valve  2  opens  inducing 
compressed  charge  impulse  through  the 
regenerator  as  the  expansion  piston 
starts  the  dawn-stroke;  valve  2  then 
cl oses. 

4.  Expansion  cylinder  exhaust  stroke: 
valve  3  opens,  regenerating  heat 
exchange  capacity  during  exhaust  gas 
expiration.  Valve  2  remains  closed; 
valve  1  is  cycling  to  support  the  next 
charge  cycle  as  expansion  piston 
reaches  top-dead-center  (TDC) . 

A  vertical  cross-section  schematic  o-f  the  expansion 
cylinder  and  piston  linkage  showing  principle  dimensions  is 
provided  in  Figure  1-2.  A  derivation  o-f  geometric 
relationships  required  to  calculate  instantaneous  expansion 
cylinder  heat  transfer  surface  area,  A^y,  volume,  V,  and 
piston  speed,  Sp ,  is  enclosed  as  Appendix  B.  Instantaneous 


expansion  cylinder  pressure,  P ,  mean  gas  temperature,  Tg,  and 


y 


volume,  V,  at  5°  crank  angle  increments  -from  the  computer¬ 
generated  thermodynamic  model  is  inclosed  as  Appendix  A, 
Table  A— 1  C13;  associated  P-V  and  Tg-V  diagrams  -for  the 
entire  cycle  are  illustrated  in  Figures  1-3  and  1-4 
respecti vely.  Other  principle  design  and  operating 
parameters  -for  the  expansion  cylinder  required  to  support  the 
heat  transfer  calculations  are  listed  in  Appendix  A. 
Appendices  A  and  B  provide  the  data  point  basis  for  the  heat 
transfer  calculations.  A  consideration  of  piston  ring  land 
area,  cylinder  sleeve  clearances,  and  clearance  volume,  Vc , 
at  TDC  is  neglected.  All  the  area  and  volumetric 
calculations  are  based  on  simple  cylinder  geometry.  Within 
the  accuracy  of  reciprocating  engine  heat  transfer 
correlations  in  general,  errors  associated  with  these 
assumptions  are  considered  to  be  insignificant. 
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FIGURE  1-1:  New  Cycle  Engine  vertical  cross-section  Eli 
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FIGURE  1-3:  New  Cycle  Engine  firing  cycle  ?-V  diagram; 
natural  aspiration  [JJ 
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CHAPTER  TWO 


FORCED  CONVECTION  HEAT  TRANSFER  IN  THE 
NEW  CYCLE  ENGINE  EXPANSION  CYLINDER 

1 .  THEORETICAL  BACKGROUND 

When  a  temperature  difference  exists  between  a  solid 
surface  and  a  moving  -fluid,  thermal  energy  is  transferred  by 
the  joint  process  of  conduction  and  convection.  Shear 
stresses  acting  on  the  fluid  will  slow  the  fluid  immediately 
adjacent  to  the  surface,  forming  a  thin  boundary  layer  of 
thickness  &.  The  heat  transferred  between  the  bulk  fluid  and 
the  body  is  transported  by  conduction  in  the  boundary  1 ayer 
as  a  function  of  fluid  movement  source  (natural  or  forced), 
flow  character  (laminar  or  turbulent),  physical  properties  of 
the  fluid,  and  body  geometry  and  roughness.  The  heat 
transfer  rate,  Q,  by  conduction  in  any  direction  x  at  any 
instant  is  determined  using  Fourier's  Law  by  a  heat 
conduction  coefficient,  K,  and  the  temperature  gradient: 

Q  =  — KA’OT  (1) 

Dx 

In  laminar  flow,  fluid  particles  flowing  parallel  to 
the  body  surface  conduct  heat  in  a  normal  direction.  At  a 
certain  transition  velocity,  turbulence  is  created  where 
particles  have  a  velocity  component  normal  to  the  bulk  stream 
direction;  the  intense  mixing  of  individual  layers  promotes 
rapid  transportaui on  of  heat  to/from  the  surface.  Even  in 
turbulent  flow,  however,  a  thin  boundary  1 ayer  exists  due  to 


ft  L1  ■  1 


-fluid  viscosity  across  which  laminar  flow  and  heat  transfer 


character 1 st 1 cs  are  maintained.  In  turbulent  systems,  the 
heat  transferred  between  a  body  surface  and  moving  fluid 
remains  limited  by  boundary  layer  thickness,  £>  ,  and  fluid 
thermal  conducti vi ty ,  K. 

An  ideal  linear  velocity  and  temperature  gradient 
across  the  boundary  layer  between  a  gas  medium  and  a  body 
surface  allows  approximation  of  the  temperature  gradient  as 
<Tg-Tw)/&;  substitution  into  equation  <1)  yields  the 
Newtonian  expression  for  steady  flow  conditions: 

q  =  -kavt 

dX 

=  -KACTq-Tw] 

& 

=  -K(Tg-Tw) 

6 

where  letting  h  x  -K  : 

& 

Q  =  hA(Tg-Tw)  (2) 

For  real  boundary  layers  of  varying  velocity  and  temperature 
gradients,  h  may  be  determined  in  some  cases  by  integration 
of  the  differential  equations  describing  tne  heat  transfer 
process.  For  more  general  and  practical  appl i cati ons , 
however,  laws  of  similarity  may  be  used  to  predict  the  value 
of  h.  According  to  the  laws  of  similarity,  it  is  sufficient 
to  determine  the  characteri stic  relations  between  independent 
dimensionless  groups  obtained,  recognizing  that  any  such 
relations  will  apply  to  all  similar  heat  transfer  processes. 
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From  the  Buckingham  pi  theorem  which  allows  reduction  of 
initial  variables  to  minimum  dimensionless  groups,  the 
•following  dimensionless  numbers  are  used  to  describe  -forced 
convection  heat  transfer  in  a  -fluid  medium  to  a  solid  body: 


1.  Nusselt  #,  Nu  =  hlc/K  (3) 

2.  Reynolds  #,  Re  =  Ul(_/V|<  (4) 

3.  Prandtl  #,  F'r  =  v../«  =  ,h.C_/K  (5) 

K  p 

where:  =  y./  p  =  momentum  diffusivity 

(kinematic  viscosity)  (6) 

«  =  K/CpP  =  thermal  diffusivity  (7) 

/i  =  dynamic  viscosity 

P  =  density 


Cp  =  specific  heat  at  constant  pressure 
lc  =  Re  characteri sti c  length 
and  in  general  form: 

Ru  =  function  of  Re,  F’r  =  f  (Re, F’r)  (8) 

In  fully  developed  turbulent  flow  systems,  analysis  of  heat 

transfer  processes  is  enhanced  by  a  simple  analogy  between  heat 

and  momentum  transfer  which  allows  prediction  of  h  from  knowledge 

of  shear  stress  (friction)  acting  on  the  fluid  adjacent  to  the 

wall.  Known  as  the  Reynolds-Colburn  Analogy,  it  is  described  by 

Rohsenow  in  the  following  manner  C2,p.lB5D: 

"the  heat  transferred  to  a  surface  divided  by  the 
maximum  amount  which  could  be  transferred  in 
bringing  the  fluid  to  the  temperature  of  the 
surface  equals  the  momentum  transferred  in 
passing  over  a  surface  divided  by  the  maximum 
amount  of  momentum  which  can  be  transferred  in 
stopping  the  fluid  relative  to  the  surface. " 

This  analogy,  exactly  true  in  parallel  flow  when  F'r  =  1  and 

reasonably  accurate  at  Fr  values  not  very  high  or  very  low, 
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may  be  used  in  either  constant  heat  -flux  or  wall  temperature 
conditions  and  may  be  expressed  mathematically  as: 


j  Colburn  -factor  =  StF 'r^'*'  =  C^/2  (9) 

where:  Stanton  #,  St  =  h/pCpU  (10) 

Cr  =  T,V(  pU^/2)  =  -friction  coe-f -f  i  ci  ent  (11) 

r  =  apparent  shear  stress  acting  on  the  -fluid 
at  the  wall 


Significant  flow  experimentation  has  determined  that  for  low 
turbulent  flow  conditions  over  a  smooth  plate  at  a  distance  x 
from  the  leading  edge: 


C,  =  . 0592 

"^flat  plate 

Similarly,  for  turbulent 
pipe  of  diameter  D: 


Re. 


flow 


conditions 


i  n 


a 


(12) 

cyl i ndr i cal 


Cf  —  • 

Tpipe  flow 

Substitution  of  these 
the  determination  of  h: 


046  ReD  ' x 
relations  into 


equation 


(13) 

(9)  allows  for 


,  flat  plate  of  length  1: 


stpr2/3  =  cf/2 

h  Pr2/3  =  .0592  (Re“,2>/2 

PCpU 

hx  =  .0296  pCpU(Rex)“*2  Pr~2/3 

=  .0296  pC  U(Re  )“>2  Pr~2/3  x  K 

P  i<  x 
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substituting  for  density,  the  equation  (6)  identity,  p  = 


.-2/3 


h„  =  .0296  K  Ux  }i£„  (Re„ )  ‘ "  Pr 

x  vk  K 

r  ecogn  l  sing  equat  i  ons  ( 4 )  and  <  5 )  l  den  t  i  ties  Re,,  =  Ux_  and 


Pr  =  >iCp/K  : 

hx  =  .0296  K  Re;<  ’  8  P r1/^ 

X 

integrating  over  length  1  to  determine  mean  : 


h.  =  1  f1  .0296  K  (U„/v,)‘8  Pr1/3  dx 

1  —  I  —  /N  rs 

1  ,  x 

J  (J 


0296  K  ( U / *v k  )  *  8  Pr  1 73 

r1  x~*2  dx 

1 

J  0 

0296  K  ( U /  v '  8  Pr  1/3 

Cx*8!1 

(J 

0.81 

037  K(Re1)*8  Pr1/3 

(14) 

A  similar  substitution  of  equation  (13)  for  turbulent  flow  in 
a  cylindrical  pipe  yields: 


hp,  pipe  flow  of  diameter  D: 


hD  =  .023  K(ReD)*8  Pr1/3‘  (15) 

D 

These  particular  relations  form  the  basis  of  many  heat 
transfer  correlations  commonly  used,  rigorously  applicable  in 
fully  developed,  steady  turbulent  flow  within  established  Re 
and  Pr  limits  where  temperature  difference,  AT,  across  the 


boundary  layer  is  not  extreme. 


2.  MODELING  CONVECTIVE  HEAT  TRANSFER  IN  RECIPROCATING  ENGINES 


In  reciprocating  engines,  the  significant  process  by 
which  termal  energy  is  lost  is  by  forced  convective  heat 
transfer  from  the  gas  medium  across  cylinder  inner  surface 
boundaries.  Large  variations  in  the  magnitude  and  direction 
of  convective  heat  transfer  exist  from  engine  to  engine  and 
point  to  point  within  an  individual  engine.  The  conditions 
influencing  heat  transfer  are  so  complex  and  locally 
different  through  an  engine  cycle  that  rigorous  mathematical 
i nterpretati on  is  invalid.  Radical  variations  in 
instantaneous  flow  patterns  allow  gas  velocities  to  be  only 
vaguely  characterized.  Non-uniformity  of  physical  properties 
in  the  gas  medium  and  surface  materials  exists.  The 
combustion  process  is  particularly  difficult  to  model  and 
further  complicated  by  compression  ignition.  Two  of  the  most 
serious  obstacles  question  the  fundamental  assumption  of 
fully  developed,  steady  turbulent  flow  and  the  concept  of  a 
finite,  instantaneous  value  of  h. 

In  operating  engines,  the  entire  cycle  is  completed 
within  milliseconds,  making  accurate  measurement  of  local 
instantaneous  heat  flux,  q,  in  itself  quite  difficult. 
Turbulent  intensity  in  four  stroke,  Cl  engine  cylinders  is 
produced  by  the  suction  velocity  of  the  induction  stroke,  the 
intensity  being  roughly  proportional  to  it.  To  optimize 
fuel-air  micro-mixing,  off-set  induction  ports  may  be 
incorporated  to  create  swirling  gas  motion.  In  the 
compression  stroke  within  a  simple  geometry  cylinder, 
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induction  turbulence  is  suppressed. 


For  swirl 


enhanced 


cylinders  with  combustion  chambers  cut  into  the  piston  crown, 
conservation  of  angular  momentum  in  the  cup  created  by 
viscous  drag  accelerates  swirl  to  significant  angular 
velocities  as  the  piston  reaches  TDC.  Radial  "squish'1 
velocities  are  created  near  TDC  and  combine  with  swirl  to 
produce  a  complicated  toroidal  motion.  In  the  expansion 
stroke,  rapid  deaccel  eration  o-f  rotational  velocities  in 
direct  injection  cylinders  are  superimposed  by  ill-defined, 
small  combustion  induced  velocities  as  hot  combustion  product 
gases  expand;  in  indirect  injection  cylinders,  prechambers 
deliberately  channel  and  exaggerate  the  charge  velocities  to 
blow  out  o-f  energy  cells  at  levels  up  to  SOOm/sec  C33.  In 
the  exhaust  stroke,  the  combustion  products  are  pushed  out  by 
the  piston.  The  accumulation  o-f  these  -factors  indicate  that 
turbulent  boundary  layers  are  not  steady,  but  are 
periodically  created  and  destroyed  in  a  manner  beyond  our 
ability  to  model.  In  addition,  because  of  the  thermal 
capacity  of  the  gases,  there  is  a  thermodynamic  phase  lag 
markedly  apparent  in  reciprocating  engines;  AT  and  q  do  not 
pass  through  zero  at  the  same  instant  indicating  values  of 
h  =  ±<»  and  highly  fluctuating  values  immediately  about  it. 

Recognizing  these  obstacles,  investigators  have 
concluded  that  within  obtainable  accuracy  it  is  still  best  to 
model  heat  transfer  as  a  quasi-steady  process  and  to 
correlate  cyclic  observations  in  terms  of  conventional  Nu , 
Re,  Pr  numbers.  Values  of  the  various  constants  and 


exponents  are  determined  in  individual  engine  tests;  a  wealth 

o-f  experimental  data  provides  guidance  tor  constant  (s) 

selection  -for  operating  engines  or  design  studies  o-f  similar 

construction.  Hohenberg  states  [4,  pg.  27881: 

"It  would  be  erroneous  to  look  -for  a  -formulation 
which  includes  all  passible  influences  on  the 
heat  transfer.  This  would  only  render  the 

calculation  more  complicated  rather  than  more 
accurate.  The  aim... must  be  to  find  out  and 
mathematically  describe  the  really  essential 
factors  influencing  heat  transfer  in  the  engine." 

If  the  time  for  change  in  bulk  stream  conditions  to  diffuse 

through  a  boundary  layer  is  small  relative  to  other 

significant  times  in  the  engine  cycle,  then  boundary  layers 

may  be  considered  quasi-steady  and  at  any  instant  be  that 

associated  with  existing  conditions.  Le  Feuvre  et  al  [53 

<2 

approximating  diffusing  time  as  &  /vk  and  modeling  gas  motion 
in  solid  body  rotation,  successfully  performed  such 
calculations  in  a  direct  injection,  swirl  enhanced,  Cl  engine 
operating  at  high  RPM  <2000  RPM)  .  As  for  the  existence  of 
phase  lag,  Annand  C6 3  states  that  the  actual  quantity  of  heat 
transferred  during  this  period  is  minute;  the  error  arising 
from  the  use  of  a  finite  h  at  each  instant  of  the  piston 
stroke  is  negligible. 

The  majority  of  the  heat  transferred  within  a 
reci procati ng  Cl  engine  is  due  to  turbulent  forced 
convection.  On  the  basis  of  experimental  data  and  theoretical 
considerations,  this  relationship  is  expressed  by  the 
formula: 

Nu  =  constant <Re) n  Pri  <16) 


Recognizing  that  Pr  .s  0.7  for  common  gas  mediums  due  to 
predominance  of  0,-,  and  ,  substitution  o-f  the  ideal  qas  law 
identity,  p  =  P/RTg,  into  an  expanded  equation  (16)  yields: 

h  =  constant  K^~nPnUnTg-nl(_r,~1  (17) 

Assuming  that  actual  gas  velocities  will  be  at  least 
approx i matel y  proportional  to  mean  piston  speed,  a  -function 
o-f  Sp  is  used  with  the  multiplicity  constant  determined  by 
engine  design.  Tg  is  used  vice  the  average,  (Tg  +  Tw)/2,  as 
the  calculation  o-f  Tg  -from  measured  instantaneous  pressure 
and  gas  law  application  already  represents  an  average  for  the 
cylinder  contents.  A  consi derati on  of  gas  di sassoci ati on  is 
neglected,  as  any  attempt  to  quantify  it  would  be  eclipsed  by 
the  previous  assumption  of  uniform  properties.  During  the 
combustion  process,  gas  product  composition  is  assumed  to 
equal  the  burned  composition.  Tables  and  Figures  of  thermal 
properties  <>».,  Cp ,  Pr ,  J'>  for  individual  gases  and  CnH2n  lean 
hydrocarbon  fuel  combustion  products  are  provided  in  Appendix 
C  C6,73.  Cyclic  local  instantaneous  surface  temperature 
profiles,  from  which  actual  measurement  of  q  is  generally 
derived,  are  demonstrated  by  Dguri  et  al  C83  to  show 
transient  magnitudes  typically  less  than  50°F  for  sooty 
surfaces  and  less  than  20°F  for  clean  surfaces.  With  the 

exception  of  exhaust  valve/port  entries  which  may  typically 
see  temperatures  greater  than  1000° F,  the  temperature  range 
exhibited  by  cylinder  sleeve,  head,  and  piston  crowns  are 
fairly  consistent  allowing  assumption  of  a  mean  Tw. 


3.  WQSCHNI'S  CORRELATION  FOR  MEAN  INSTANTANEOUS  HEAT 
TRANSFER  COEFFICIENT,  hg,-  IN  Cl  ENGINES 

A  recognized  industry  standard  tor  -four  stroke,  direct 
injection,  quiescent,  Cl  engine  heat  transfer  calculation  is 
Woschni's  correlation,  stated  in  U.S.  units  as  follows  C91: 

he  =  15.480  P*8CC1Sp+C2<VTgc/F'cVc)  (P-Po)  3  •  8B~*  2Tg“‘ 53  (18) 
where  he  =  BTU/hr-f t2-°R 

F'c ,  Vc,  Tgc  =  i  nstantaneous  values  at  ignition  start 
Po  =  motoring  pressure 

Based  on  suppressed  swirl  engine  designs,  the  values  of 
constants  and  C2  ares 


1. 

induction  stroke  C,  =  6.  18 

C-,  =  0 

(18a) 

n 

S-  m 

compression  stroke  Cj  =  2.28 

Co  =  0 

dL 

( 18b ) 

•— 1  m 

expansion  stroke  (motoring)  Cj 

—  O  '"■O 

—  A—  •  a—  O 

( 18c) 

Co 

Am 

=  0 

4. 

combustion-expansion  stroke  Cj 

—  o  '"JO 

“ “  Amt  ■  dm  O 

( 18d) 

C2  =  5.9054  x  10“ J  f t/sec-°R 

5. 

exhaust  stroke  Cj  =  6.  18 

Co  =  0 

A- 

( 18e) 

From  equation  (17),  a  n  exponent  of  0.8  was  determined  by 
Woschni  by  successively  changing  one  parameter  at  a  time  over 
maximum  possible  range.  A  consensus  opinion  of  all  recent 
published  correlations  is  that  0.7  6  n  «=  0.8,  reflecting 
agreement  with  the  laws  of  similarity.  Woschni  developed  the 


motoring/f iring  correlation  and  scavenging  correlation  on 
suppressed  swirl  design  cylinders  stating  L9,p.3071D: 


"The  gas  motion  generated  during  the  induction 
phase  dies  away  with  time,  owing  to  internal 
-friction.  In  the  compression  stroke,  however, 
dependent  upon  the  piston  or  cylinder  cover 
shape,  an  additional  gas  motion  can  be 
induced  C squi sh-swi r 1 3  whose  intensity  is  also 
proportional  to  piston  speed.  In  this 
particular  case,  the  engine  involved  had 
direct  injection  and  straight  -forward 
combustion  chamber  walls,  so  that  decay  due  to 
internal  -friction  prevailed." 

A  Tg  exponent  is  due  to  Woschni  's  assumption  o-f  air  and 
gas  product  ,u.(Tg)  «  Tg‘6and  K(Tg)  c<  Tg*7^  at  n  exponent  =  0.8 
(see  equation  (17)).  To  approximate  lll-de-fined  velocity 
increase  due  to  combustion,  the  difference  P-Po  which,  unlike 
Tg-To,  sees  a  rapid  decay  upon  combustion  completion,  is 
used,  multiplied  by  instantaneous  cylinder  volume,  V,  and 
related  to  the  charge  mass  at  beginning  of  combustion, 

mc  =  pc'7c/RT^c* 

It  is  at  this  point  that  Woschni  makes  a  critical 
simplifying  assumption.  Recognizing  the  significant  effect 
of  luminous  flame  radiation  heat  transfer  in  Cl  engines, 
Woschni  arbitrarily  accounts  for  combustion  radiant  heat 
transfer  by  affixing  the  value  Cj  to  equate  actual  measured  Q 
integrated  over  the  stroke  with  Q  combustion  determined  by 
the  correlation.  Woschni  states  C9,p. 30753: 


"There  is  no  doubt  that  such  an  influence  occurs 
in  the  engine,  and  it  is  taken  into 
consideration  in  lumped  form  in  the  combustion 
term  as  defined  here,  but  it  seems  impossible  to 
the  author  to  separate  it  from  the  convective 
term.  The  only  question  that  remains  to  be 
considered,  therefore,  is  the  error  made  regard 
to  the  tendency ...  for  the  radiation." 


In  essence,  for  application  to  the  New  Cycle  Engine 
expansion  cylinder  operating  at  significantly  lower  peak 
pressure  conditions,  laws  of  similarity  allow  significant 
extrapol ati on  of  Woschni's  correlation  for  forced  convection 
heat  transfer  only.  The  correlation  must  be  corrected  to 
account  for  the  fact  that  radiant  heat  transfer  is  not  a 
function  of  <P— Po)‘®. 

4.  CALCULATION  OF  ORDER  OF  r-IAGNITUDE  REDUCTION  IN 

COMBUSTION  FORCED  CONVECTION  HEAT  FLUX 

The  representative  firing  cycle  P-Tg  versus  crank  angle 
plot  published  by  Woschni  for  natural  aspirating  operation  is 
shown  as  Figure  2-1  C93.  The  range  of  additional  principle 
operating  parameters  over  which  the  motor i ng/f 1 r i nq 
correlations  were  developed  wel 1  representive  the  essential 
parameters  of  New  Cycle  Engine  (NCE)  expansion  cylinder 
operation  C91: 

Parameter  Woshni  (W)  engine  NCE  Expansion  Cylinder 

RPM  900  max  SI 8. 3 

5p  9  m/s  max  C29.6  f/sl  33.33  ft/sec 

« 

Supercharge  3  max  1.0 

pressure 
ratio 

Tw  nominal  500° K  C9009R1  900® R 

L,B  unknown  1.22204  ft 

Bore  and  stroke  dimensions  were  not  provided,  however  from 
known  maximum  data  point  values  Sp  =  9  m/sec  and  NRpM  =  900, 
the  approximate  stroke  length  L  may  be  determined: 


■>£> 


L  «  60  Sp  sj  (60)  (9)  X  0.3  m  CO.  98  f  t  ] 

2Np,pM  (2)  (900)  approx. 

For  conventional  Cl  engine  geometry,  bore  diameter  B  is 

approx  i  matel  y  equal  to  L;  there-fore  -for  effective  comparison 

—  _  'T' 

of  be  proportional  to  B  " ^  in  forced  convection  heat 
transfer,  consi derat i on  of  bore  diameters  of  similar 
magnitude  dimensions  may  be  neglected.  The  significant 
difference  between  combustion-expansion  stroke  profiles  ares 

1.  An  order  of  15  magnitude  increase  in  peak 
combustion  pressure  in  Woschni's 
conventional  high  rc ,  Cl  engine. 

2.  A  reduced  Tg  decay  rate,  thus  increased 
<Tg-Tw)  farced  convective  heat  transfer 
potential  in  the  New  Cycle  Engine 
expansion  cylinder. 

3.  A  delayed  fuel  injection/ignition  start 
in  the  New  Cycle  Engine  expansion 
cyl i nder . 

Assuming  that  the  magnitude  of  mean  instantaneous  gas 
velocities  through  the  expansion  strokes  are  roughly  equal, 
the  approximate  ratio  of  forced  convective  heat  transfer  loss 
in  the  respective  engines  may  by  obtained  by  determining  the 
mean  ’ratio  of  instantaneous  heat  flux,  q,  over  a  90®  crank 
angle  domain  from  ignition  start.  A  selection  of  90°  is 
based  on: 

1.  Establishing  a  normalized  basis  in  which  to 
compare  the  order  of  magnitude  reduction  in 


forced 


convection  versus 


radi ant 


heat 


transfer;  the  crank  angle  domain  over  wnich 
luminous  flame  radiation  is  significant  in 
medium  speed  Cl  engines  less  than  100U  ft PH 
is  nominally  90°  from  combust  1  on  start  181. 

2.  At  90°  the  slope  of  pressure  and 
temperature  decay  transients  induced  by 
combustion  in  both  engines  become  small. 

Setting  the  radiant  heat  transfer  constant  C2  =  0  in  equation 
(18d)  : 

qW/qNCE  w  ^e (Tg— Tw) w/he (Tg-Tw) NCE 


qW/qNCE  c< 


CP*8Sp*8/B*2Tg*533CTg~Tw3w 
CP* 8§p*  S/B“ 2Tg*  53 3  CTg-Tw3NrE 


Neglecting  B  consideration  and  assuming  equal  magnitude  gas 
vel  oci t i es: 


qw  CP*8/Tg*003CTg-Tw3w 

_  «  - - - — -  (19) 

qNCE  CP*8/Tg*'J'iCTg-Tw3NCE 


where  7w  =  900° ft 


Table  2-1  summarises  the  calculations  based  on  Figure  2-1  and 
Appendix  I  data  at  5°  increments  from  combustion  start.  The 
mean  heat  flux  reduction  ratio  may  then  be  determined  by 
numerical  integration  using  Simpson's  rule; 
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TABLE 


1 


Calculation  of  order  of  magnitude  reduction  in  combustion 
forced  convection  q  from  Woshni  s  test  engine  to  New  Cycle 
Engine  expansion  cylinder  operation  at  5°  crank  angle 
1 ncrements 


Crank  anglt 

Noshni 

tngim  (M) 

NCE  Expansion  Cylinder 

V%e 

past  ignition 

Tg 

P 

Tg 

P 

Tg 

aquation  (19) 

start 

kp/cr 

•K 

at* 

•R 

at* 

•R 

0 

40 

873 

38.71 

1373.0 

4.13 

1896.6 

4.46 

5 

60 

1400 

58.07 

2320.0 

4.23 

2252.6 

9.14 

10 

63 

1323 

60.97 

2743.0 

4.09 

2371.2 

9.29 

IS 

60 

1600 

38.07 

2880.0 

3.92 

2864.1 

8.68 

20 

53 

1650 

53.24 

2970.0 

3.76 

3136.0 

7.94 

23 

47 

1630 

43.49 

2934.0 

3.26 

3064.2 

• 

7.92 

30 

40 

1600 

38.71 

*880.0 

2.89 

2998.7 

7.68 

33 

33 

1540 

31.94 

2772.0 

2.39 

2939.1 

6.66 

40 

26 

1480 

25.16 

2664.0 

2.34 

2884.8 

6.2 

43 

22.3 

1440 

21.18 

2392.0 

2.13 

2835.2 

5.76 

30 

20 

1400 

19.36 

2520.0 

1.96 

2790.0 

5.63 

35 

17 

1363 

16.43 

2430.0 

1.81 

2748.8 

5.16 

60 

13 

1330 

14.32 

2394.0 

1.69 

2711.3 

4.92 

63 

12.3 

1305 

12.1 

2349.0 

1.38 

2677.1 

4.46 

70 

11 

1280 

10.63 

2304.0 

1.49 

2646.0 

4.17 

75 

10 

1233 

9.68 

2259.0 

1.42 

2617.8 

3.97 

80 

9 

1230 

8.71 

2214.0 

1.33 

2392.3 

3.73 

85 

8.5 

1213 

8.23 

2187.0 

1.3 

2369.2 

3.68 

90 

3 

1200 

7.74 

2160.0 

1.23 

2348.4 

3.59 

40 


u  • 


E  S  TIMA T ION  OF  MEAN  GAS  VELOCITY  MAGNITUDE,  U 


5. 1  CONSTRUCTION  OF  CHARGE  IMPULSE  VEUOC I T V , U 1 , T PANS I EN T 

Acceptance  of  the  general  applicability  o-f  Woschni's 
correlation  to  the  New  Cycle  Engine  expansion  cylinder 
balances  on  the  assumption  o-f  comparable  gas  velocity 
magnitudes.  While  geometric  -flow  patterns  are  indeterminate 
in  preliminary  design,  the  magnitude  o-f  the  charge  impulse 
into  the  expansion  cylinder  leaving  the  regenerator  port  can 
be  ascertained.  Recognising  the  inevitable  deaccel er at  1  on  o-f 
this  impulse  in  a  simple  geometry  cylinder  on  the  expansion 
stroke,  a  conservative  estimation  of  the  mean  bulk  gas 
velocity  at  the  walls,  Uexp,  through  the  expansion  stroke  can 
be  made  with  minimum  error.  Substituting  Uexp  directly  into 
the  steady,  fully  developed  turbulent  flow  equations  (14)  and 
(15)  for  flat  plate  and  pipe  flow  should,  if  the  magnitudes 

of  cylinder  velocities  are  in  fact  similar,  be  reasonably 

• 

close  to  he  calculated  from  Woschni's  correlation.  To  reduce 
uncertainties,  the  relationships  will  be  first  compared  in 
motoring  conditions  from  TDC  to  45°,  prior  to  ignition  start. 

Dent  et  al  C1QD  successfully  used  a  similar  method  in 
approx i mati ng  forced  convective  heat  transfer  loss  in  a  four 
stroke,  direct  injection,  high  swirl,  Cl  engine  operating  at 
peak  TDC  swirl  ratios  equal  to  20.  Predictions  using 
Woschni's  motoring  correlations  failed  uniformly  low  compared 
to  measured,  spatially  averaged,  instantaneous  heat  flux. 
This  is  expected,  considering  the  central  assumption  of 
suppressed  swirl  inherent  in  application  of  Woschni's 


correl ati on 


By  substituting  directly  into  the  + 1  at  plate  h 


equation  (14),  measured  instantaneous  angular  gas  velocity  uy 
at  lc  equal  to  bore  c  1  rcumf erence  (2trr>,  remarkable  agreement 
was  obtained.  Agreement  was  not  as  good  at  tiring  conditions 
due  to  the  transient  combustion  processes. 

From  ideal  gas  law,  at  a  given  angle  A  in  the  expansion 
stroke,  charge  mass,  m,  in  the  expansion  cylinder: 

F'V  (2116.224  lbf /f t2-atm)PV 

m6  lbm  ~  T  -  I 

RTg  (53. 90) Tg 


Taking  5° 
RPM;  the 


=  39.262  PV 
Tg 

crank  angle  as  the 
charge  mass  rate, 


(21 ) 

unit  increment  of  time  at  818.3 
entering  the  expansion  cylinder 


is: 


m6  lbm/sec 


Am(5° ) 
At  <5°  ) 


Am  <5°  ) 

1 _ 5 _ 

(NRpM/60)360® 


=  Am (5° )  (22) 

.0010184  sec 

The  continuity  equation  for  impulse  flow  leaving  the 
regenerator  port  cross-section  area,  Ap  ,  equal  to  0.80  ft*' 
is: 

=  pl/iAp  (23) 

Therefore,  by  rearranging  equation  (23),  impulse  velocity 
leaving  the  regenerator  port,  Ui ,  can  be  determined: 
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Ui 


m9RT  g 


m^ (53. 90) Tg 


ft/sec 


PA, 


(2116.224  lbf /f t^-atm) P ( . 080) 


.31337  Tg 
P 


(24) 


A  summary  of  the  calculations  using  instantaneous  F‘,  Tg,  V 
values  -from  Appendix  A  at  5°  increments  is  provided  in  Table 
2-2.  Charge  mass  as  a  function  o-f  expansion  stroke  crank 
angle  is  plotted  in  Figure  2-2,  indicating  that  the  main 
impulse  from  the  compression  cylinder  via  the  regenerator 
exit  port  completes  10°  past  ignition  start. 

At  45°  combustion  start  "/,mtotal  ~  2 •  10~t!  =  39.37. 

2 . 4537 x 1 0"2 


At  55°  main  impulse  completion 

7.m 


total 


28894x  10  *-  =  93.37. 


2.4537x10  “■ 


Ui  as  a  function  of  expansion  stroke  crank  angle  is  plotted 
in  Figure  2-3,  indicating  a  peak  Ui  equal  to  519.95  ft/sec 
near  fuel  injection  start  <30°). 
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TABLE  2-2 


Calculation  of  charge  Ui  leaving  the  regenerator  port  into 
the  New  Cycle  Engine  expansion  cylinder  at  5°  crank  angle 
1 ncrements 


Angle 

m 

m 

Ui 

6° 

1  bm 

A1  bmE— 3 

1 bm/sec 

■f  t  /sec 

eq. (21 ) 

eq. (22) 

eq. (24) 

TDC 

0 

0 

0 

0 

5 

1 . 5550E-4 

0 . 1 555 

0. 1527 

41.82 

10 

7. 7090E-4 

0.6154 

0. 6043 

133.  14 

15 

2. 1309E-3 

1 . 360 

1 . 3275 

236.92 

20 

4. 5821E-3 

2.451 

2. 4069 

352.69 

25 

8.3113E-3 

3.  730 

3 . 66 1 8 

458. 13 

30 

1 . 2936E-2 

4.625 

4.5413 

519.95 

35 

1 . 7373E-2 

4 . 437 

4.357 

499.  15 

40 

2. 049 IE— 2 

3.  118 

3 . 06 1 4 

382. 68 

45 

2. 1905E-2 

1 . 414 

1 . 3882 

202. 03 

50 

2.2887E-2 

0 .9816 

0.9641 

162.73 

55 

2.2889E-2 

0. 0028 

0. 00275 

0 . 55 

60 

2 

0. 0426 

0. 04183 

9.  72 

65 

2. 2992E-2 

0. 0596 

0.05852 

15.54 

70 

2. 3056E-2 

0 . 064 1 

0. 06294 

18.83 

75 

2 . 3307E-2 

0.251 0 

0 . 2465 

81 . 44 

80 

2.3512E-2 

0. 2049 

0. 2012 

72.  76 

85 

2. 3680E-2 

0 . 1 687 

0. 1657 

65.12 

90 

2.3820E-2 

0 . 1 397 

0. 1372 

58.15 

95 

2. 3937E-2 

0. 1167 

0.1146 

52.  03 

100 

2. 4035E-2 

0. 0981 

0. 09633 

46.  56 

105 

2. 41 18E-2 

0.0827 

0.08121 

41.54 

no 

2.4188E-2 

0.0702 

0. 06893 

37.  12 

115 

2.4247E-2 

0.0596 

0 . 05852 

33.  0 

120 

2. 4298E-2 

0.0511 

0. 05018 

29.84 

125 

2.4342E-2 

0.0436 

0.04281 

26.  13 

130 

2. 4379E-2 

0.0372 

0 . 03653 

23 . 04 

135 

2.44UE-2 

0. 0320 

0.03142 

20.41 

140 

2.4438E-2 

0.0270 

0. 02651 

17.67 

145 

2. 4462E-2 

0.0239 

0. 02347 

15.99 

150 

2.4480E-2 

0 . 0 1 83 

0.01797 

12.48 

155 

2. 4496E-2 

0. 0158 

0.01552 

10.  94 

1 60 

2. 4509E-2 

0.01 30 

0.01276 

9.  12 

165 

2. 4519E-2 

0 . 0 101 

0.00992 

7.  15 

170 

2. 4527E-2 

0.0074 

0. 00727 

5.28 

175 

2. 4532E-2 

0 . 005 1 

0. 00501 

3.  o5 

BDC 

2.4537E-2 

0.0052 

0. 0051 1 

3.  70 

45 


p 


5.2  VERIFICATION  OF  TURBULENT  FLOW 

To  verify  turbulent  flow  conditions  leaving  the 

regenerator  port,  calculated  Re  at  30°  using  lc  equal  to  the 

regenerator  port  cross-section  minor  dimension  (depth)  of 
0.75  inch  is: 

Uilc  Ui 1 c  •  Ui 1 c  P 

Re-,00  =  =  -  -  - -  (25) 

vk  }L/p  n  RTg 

From  Appendix  C,  Table  C-l ,  >i<Tg)  and  Cp(T)  for  air  is: 

Ibm/ft-sec  =  2. 163x 10~7 (Tg) * 64b  (26) 

(Cl) 

C  BTU/1 bm-hr  =  0.224  +  0.262xl0~4  Tg  (27) 

(C2> 

Substituting  instantaneous  values  of  F',  Tg  from  Appendix  A 
and  )i(Tg)  identity  equation  (26)  into  equation  (25): 

P7.-.0  -  6.10171  atm 
2194. 29°R 

(519.95) (0.75/12) (2116.224  lbf/ft2-atm) (6.10171) 
(2.163xl0-7) (2194.29) ' 645  (53.90) (2194.29) 

1.1476x10  turbulent 


Tg30°  ~ 


Re30° 
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5.3 


The  calculation  of  the  mean  impulse  velocity,  Di  ,  -for 
the  purpose  of  Uexp  approximation,  is  correctly  restricted  to 
the  main  impulse  domain  TDC  to  55°  only,  neglecting  residual 
10%  charge  mass  regenerator  bleed.  Combustion  induced 
perturbations  during  crank  angles  45°-55°  may  be  neglected. 
Ui  may  then  be  calculated  using  numerical  integration  by 
trapezoidal  rule  and  is  equal  to: 


‘ft/sec 
TDC— 55° 


_1_  f‘ 

£55°  J i 


( £5°  /^i)C  U  i  -j-qq  ^SUi^o  ■^2Ui^Qo...+2\Jig^o"*‘^i|j52o*l 


=  <£5°/2) C5979. 033  =  271.77  ft/sec 

£55° 

A  reasonable  assumption  must  be  made  as  to  the  expected 
value  of  Uexp  over  the  entire  stroke.  The  principle  factors 
causing  charge  velocity  deaccel er ati on  include  cross-secti on 
area  expansion  exiting  the  regenerator  port,  internal 
friction,  cylinder  volumetric  expansion,  and  angular  momentum 
deaccel erati on  of  any  existing  rotating  flow  patterns. 
References  C5,  10,  11,  12,  133  unanimously  agree  that  in  swirl 
enhanced  engines,  peak  instantaneous  swirl  at  TDC  is  much 
greater  than  the  mean  swirl  valves  over  the  stroke. 

Le  Feuvre  et  al  C53  calculated  for  a  4.5  inch  bore  direct 
injection,  Cl  engine  operating  at  2000  RPM,  that  10%  velocity 


4<? 


decrease  occurred  each  stroke  due  to  viscous  dissipation 
alone.  A  conservative  estimate  of  Uexp  =  ( Oi jDC-550 ) is 
therefore  assumed;  selection  o-f  this  value  should  reduce  any 
error  in  either  direction.  It  is  noted  that  UiTDC_BDC  -from 
Table  2-2  in  firing  conditions  equals  102.9  ft/sec. 
Combustion  induced  velocities  due  to  charge  gas  expansion  in 
spark  ignition  <SI)  engines  have  been  calculated  by  Annand 
C33  to  be  on  the  order  of  33  ft/sec,  with  expected  velocities 
in  Cl  engines  less  due  to  diffusion  flame  conf i gurat i on . 

Thus  ,  1  et  Uex p  =  ( Ui  -pQQ _ 55°  ^  /2  =  1  --'5 .9  ft/ sec  ■ 

6.  COMPARISON  OF  hg  CALCULATED  FROM  WOSHNI'S  CORRELATION 
AND  THEORETICAL  MODELS 
6. 1  EXPANSION  STROKE;  MOTOR I NS  CONDITION 

Using  the  above  first  approximation  of  Uexp,  estimates 
of  he  from  theoretical  equations  (14)  and  (15)  for  steady, 
fully  developed  flow  may  be  reached.  Expanding  pipe  flow 
equation  (15)  at  lc  =  B  and  U  =  Uexp  yields: 

hD  =  .023(K/B)ReD' 8Pr 1/3  =  . 023 (K/B) CUexpB/ >k 3  * 8Pr 1 /3 

Substituting  the  identities  >k  =  pj  p  ,  K  =  Cp/R/Pr,  and 
P  =  P/RTg: 

hD  =  .  Q23Cp/iCPUexp  / ^iRlg  1 *  8B- *  2Pr-2/3 

Substituting  the  Cp(T)  and  >i<Tg)  identities  from  equations 
(27)  and  (28)  and  appropriate  conversions: 
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hD  =  . 023C.224+.262x 10  4Tg 3 C2. 163x 10~7 (Tg ) ’ 645 (3600sec/hr ) 1* 

C (2116.224  1  bf  /-f  t2-atm  P  Oexp  )  /  (2.  163x  10_7Tg*  645RTg  )  1  *  8* 


_  O  _  O  /  T 

B  ' iPr  ^/- 


which  simplifies  at  known  R,  Uexp,  Pr  values  to: 


BTU/hr— 


=  9.194112.24+2.62x10  4Tgl- 


(135.9  P ) • 3 


B'  ^  T g  ‘ 


An  identical  expansion  of  the  -flat  plate  equation  (14)  at 


lc  =  irB,  expansion  cylinder  ci rcum-f erence ,  yields: 


BTU/hr-f  t"--°R 


=  11. 674C2. 24+2. 62x 10  4Tgl- 


(135.9  P ) • 1 3 


B ’ 2  Tg‘671 


A  summary  o-f  h  calculations  for  the  charge  induction 
monitoring  angles,  TDC-450 ,  comparing  equations  (29)  and  (30) 
to  Woshni's  motoring  correlation  equation  (18c),  is  provided 
in  Table  2-3.  The  agreement  is  quite  good,  even  with 
acceptance  of  reasonable  error  in  Dexp  approximation. 


TABLE  2 


Comparison  of  Fie  calculated  by  Woshni's  correlation  and 
theoretical  models  during  charge  induction  into  the  New  Cycle 
Engine  expansion  cylinder 


Angle 

Fe 

hD 

hl 

8° 

eq . (18c) 

eq. (29) 

eq . ( 30 

5 

16.994 

15.221 

19.477 

10 

20. 191 

18.143 

n-r  i  i  cr 

O  a  Li  J 

15 

24. 006 

2 1 . 593 

27.631 

20 

28.209 

25.396 

32.496 

25 

32 . 029 

28.839 

36.902 

30 

34 . 236 

30 .793 

39. 403 

35 

33 • 888 

30 .414 

38.918 

40 

31 . 168 

27.897 

35.697 

45 

27. 164 

24. 251 

31. 037 

h  in  BTU/hr — f t2-0 R 


6.  2 


EXHAUST  STROKE 


Attempts  to  directly  estimate  the  mean  bulk  gas 
velocity  at  the  walls  during  the  exhaust  stroke,  Oexh,  in 
light  of  preliminary  design  restrictions  and  post-combustion 
ambiguities,  are  not  valid. 

An  indirect  method  allowing  h  calculation  is  to  correct 
equation  (29)  by  a  -factor  equivalent  to  the 
scavengi ng/motor i ng  Oexh  correction  observed  by  Woschni  -for 
suppressed  swirl  geometry,  direct  injection,  Cl  engine 
cylinders.  From  equations  18a  and  18c: 

Oexh  correction  factor  =  (6.18  Sp/2.28  5p)-8  =  2.22 

When  the  regenerator  exit  valve  3  opens,  the 

significant  flow  pattern  will  be  parallel  to  the  bare  axis, 
up  the  sleeve,  and  through  the  regenerator  port;  flow 
direction  is  therefore  analogous  to  pipe  flow.  From  Appendix 
C,  the  identities  >i(Tg)  ,  Cp(Tg),  Pr  (Y)  ,  and  J'(Tg)  for  lean 
CRH2n  hydrocarbon  fuel  combustion  products  at  P  =  1  atm, 
i  =0.75  may  be  calculated: 

1 . 47  x 10~ 7 ( Tg ) *7  1 . 47x 10-7 (Tg ) '7 

>l  lbm/ft-sec  =  -  =  - 

1  +  .027#  1  +  .027 (.75) 


=  1 . 441 x 10~7 (Tg ) *  7  (31) 

(C3) 

C  BTU/lbm-oR  =  0.247  +  0.243xl0-4  Tg  (32) 

P  (C5> 

Pr(J')  =  0.05  +  4 .  2  ( if  —  1 )  -  6.  7(3-1)  2  (33) 

(C4) 

J'(Tg)  z  Appendix  C,  Figure  C-l  (34) 


r 


Applying  the  Uexh  correction  factor  =  2.22  and  the  identities 
^.(Tg)  ,  Cp<Tg)  defined  above  to  pipe  -flow  equation  (29) 
yields: 

hD  =  2. 22 ( . 023)  ( -  247+. 243x 10~4Tg )  ( 1 . 441 x 10~7Tg *  7 )  (3600sec/hr ) * 
C  (2116.224  lbf/ft2-atm  P  Oexp ) / < 1 . 44 1 x 10~7Tg * 7RTg > J ’ B* 

B  *  “Pr  't/'> 


which  simplifies  at  known  R‘cp,  Uexp  values  to: 


=  14.835:2.47  +  2.43x10  HTg] 


BTU/hr— f  t^-°  R 


135.9  P 
B-2Tq-66 


Pr-2/3 


From  Appendix  A;  P  =  1.05  atm  and  2344. 58°K  s=  Tg  ±  2439. 06°R 

TDC  BDC 

through  the  exhaust  stroke.  A  comparison  of  h  calculated 
from  equation  <351  and  Woschni ‘s  exhaust  stroke  correlation 
equation  (18e)  follows: 


crank  angle  Tg  J'(Tg)  F’r(J')  h0  hQ 

°R(°K)  eq. (34)  eq. (33)  eq. (18e)  eq. (35) 

BDC  2439.06  1.275  0.698  17.587  17.062 

( 1 355 ) 

TDC  2344.58  1.280  0.701  17.959  17.331 

(1303) 

h  in  BTU/hr-f t2-0 R 

The  above  series  of  calculations  suggest  that  to  a  first 
approximation,  application  of  Woschni  's  correlation  in 


determining  forced  convective  heat  transfer  loss  for  the  New 
Cycle  Engine  expansion  is  valid. 


CHAPTER  THREE 


CALCULATION  OP  UNCQRRECTED  HEAT  TRANSFER 
LQSS/HEAT  RELEASE  PER  CYCLE,  Tt 

Table  3.1  summarises  the  calculation  o-f  total 

instantaneous  heat  transfer  rate,  Q,  through  a  complete  36 0° 

firing  cycle  of  the  New  Cycle  Engine  expansion  cylinder  using 

Woschni's  correlation  for  h  .  For  calculation  of  h„  in 

e  e 

combustion,  motoring  p  is  generated  graphically  in  Appendix 
A,  Figure  A- Is 

QgTU/sec  =  ^e^HT (Tg-Tw) ( 1  hr/3600  sec)  (36) 

Due  to  the  establishment  in  preliminary  design  of  a  simple 
cylinder  geometry,  AHT  as  a  function  of  crank  angle  may  be 
determined  directly  from  instantaneous  V  values;  defining 
i nstantaneous  gas  volume  height  as  L^  and  neglecting  cylinder 
sleeve  clearance: 

v  ■  L« 

Art  =  2rrB2  +  TT&Le  =  2itB2  +  irDC  4V/:rB2  J 
4  4 

which  at  B  =  1.22204  ft: 

Aht  ^  =  2.346  +  3.273CV3  (37) 

ft* 

_  * 

The  distribution  of  h0  and  Q  as  a  function  of  crank  angle  is 
plotted  in  Figure  3-1.  It  is  immediately  clear  studying  the 
he  distribution  over  the  combustion  stroke  55®  to  BDC,  that 
Woschni's  correlation  fails  to  correctly  extrapolate  luminous 


flame  radiant  heat  transfer  loss  to  low  rc  operating 
conditions.  No  increase  in  h0  occurs  due  to  the  predominant 
influence  of  increased  temperature:  associated  gas  layer  p 

decrease  and  ft.  increase  overcome  the  effect  of  K  increase, 
thus  lowering  he,  modeled  by  laws  of  similarity  solely  on 
farced  convective  heat  transfer.  The  low  rc  and  delayed  fuel 
injection  operating  conditions  effectively  eliminate  any 
remarkable  pressure  rise  in  the  cylinder,  already  at  45° 
down-stroke  at  ignition  start.  Recalling  Woschni's  coupling 

O 

of  radiant  heat  transfer  in  combustion  to  C-,<F— F'  )  ,  the 

jl  (j 

desired  "lumped  form"  term  effectively  disappears. 

Uncorrected  total  heat  transfer  loss  per  cycle,  Q-j- ,  is 
determined  by  numerical  integration  under  the  Q  curve  using 
Simpson's  rule: 


L 


BTU/cycle 


~  ^01850+  ^190°  + 


1 1 7  0 c 


1 7! 


BDC' 


<38 ) 


=  -001018417141.6133  =  2.42434  BTU/cycle 


From  Table  2-2,  total  charge  mass  per  cycle  equals  .024537 
lbm.  At  known  specific  heat  release  per  cycle,  A  heat,  equal 
to  500  BTU/lbm  charge,  the  ratio  of  total  heat  transfer 
loss/heat  release  per  cycle,  Tly,  equals: 


a. 


2.42434  BTU 


uncorrected  (500BTU/ 1 bm) (. 024537  lbm/cycle) 


=  19.7617. 
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A 


The  inability  of  Woshni  s  correlation  to  account  -for  radiant 


heat  transfer  in  low  rc  operation  is  clearly  demonstrated  by 

repeating  the  numerical  integration  with  the  combustion  term 

in  equation  (18d)  eliminated  (Co  =  0): 

Simpson's  Multiplier  =  7038.088 

Qt  =  2.3892  BTU/cyc 1 e 

1lT  =  19.474% 

uncorrected 

Considering  the  accuracy  of  the  correlation  in  general ,  for 
preliminary  design  prediction  purposes  setting  =  0  is 

considered  valid. 

To  correct  for  luminous  flame  radiation  heat  transfer 
loss,  a  value  of  A^^NCE  radiation  must  be  determined  which, 
in  comparison  of  clw/c,NCE  * Q|i~ced  convection  calculated  in 
Chapter  Two,  will  allow  first  appro;; i mat i on  of  the  required 
tlj  correction  factor.  As  there  exists  no  theoretical  basis 
for  coupling  magnitude  of  luminous  flame  radiation  to 

Q 

<P-P0) *  ,  an  attempt  to  determine  the  desired  ratio  from 

s 

<P-PD> '  profiles  must  be  disregarded  in  favor  of  calculation 
founded  on  known  principles  of  radiation  heat  transfer  in 
diffusion  flame  combustion. 


TABLE  3-1 


Calculation  o-f 

• 

cycle  Q  using 

i ncrements 

crank  angle 

0 

e 

charge  induction 

New  Cycle  Engine  expansion  cylinder  tiring 

Woschni's  correlation  at  5°  crank  angle 

?g-Tw  Q 

BTU  ft2  °R  BTU 

2  SSC 
hr— f  t  — °R 

eq<18)  eq(37)  Tw=900°R  eq<36> 

equation  ( ISc) 

5 

16.944 

2.357 

1 249 . 04 

13.856 

10 

20. 191 

2 . 390 

1259.25 

16.878 

18 

24. 006 

2.445 

1295.73 

21.128 

20 

28.209 

m  JAA 

1 325 . 26 

26. 190 

25 

32.029 

2.618 

1 329 . 7 1 

30.972 

30 

34 . 236 

2.734 

1294.29 

33. 652 

35 

33. 888 

2.867 

1216.61 

32.834 

40 

31.1 68 

3.  017 

1 1 10.57 

29. 009 

45 

27. 164 

3.  181 

996.6 

23 . 924 

combustion-expansion 

equation  < lSd ) 

50 

26.086 

3 .  358 

1351.6 

32. 888 

55 

23.969 

3.548 

1671 . 21 

39. 445 

60 

22.066 

3.  741 

1964. 13 

45 • 038 

65 

20.853 

3.  944 

2236. 02 

51.083 

70 

18.710 

4.152 

2164. 15 

46. 700 
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TABLE 


1  (CQNT ' D) 


:  angle 

h- 

I 

<r 

Tg-Tw 

Q 

0 

6 

BTU 

-ft*- 

°R 

BTU 

sec 

hr-+ti-°R 

75 

17. 084 

4.  362 

2098.73 

43.444 

80 

15. 742 

4.573 

2039 . 1 1 

40.775 

85 

14.632 

4.783 

1984.75 

38. 584 

90 

13. 715 

4 . 990 

1935. 18 

36. 789 

95 

12.841 

5.  192 

1889.98 

35.002 

100 

12. 099 

5.  388 

1848. 18 

33.473 

105 

11.457 

5.577 

1811.25 

32. 148 

110 

10.865 

5.757 

1 777 . 09 

30.877 

115 

10.438 

5.  927 

1746. 03 

30. 006 

% 

120 

10.015 

6.  087 

1717.83 

29. 089 

125 

9.652 

6 . 236 

1692.29 

28.361 

130 

9 . 329 

6) » -^7  o 

1669.22 

27.567 

135 

9.052 

6.498 

1648.44 

26. 934 

140 

8 . 830 

6.611 

1629.82 

26. 428 

145 

8.  632 

6.711 

1613.21 

25.959 

150 

8.  467 

6.797 

1598. 51 

25.554 

155 

8.334 

6.871 

1585.62 

25. 221 

160 

8.229 

6. 930 

1574. 46 

24.941 

165 

8.  151 

6.973 

1564.96 

24.725 

170 

8.086 

7.011 

1557.07 

24.517 

175 

8.044 

7 . 03 1 

1550.74 

24 . 363 

B  DC 

8.028 

7 . 037 

1545.46 

24 . 260 
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TABLE 


1  <CQNT ' D  > 


crank  angle 

Re 

aht 

7g-Tw 

Q 

0 

e 

BTU 

f  t*- 

°R 

BTU 

sec 

hr— ft2-°R 

exhaust 

equation  (18e) 

BDC 

17.587 

7 . 037 

1539.06 

52.909 

185 

1 7 . 604 

7.031 

1534.63 

52.763 

190 

17.619 

7.011 

1530. 70 

c  *“>  cr  -t 
J4.  •  U* 

195 

17.634 

6.  978 

1526.81 

52. 187 

200 

17.649 

6.  930 

1522.95 

51.741 

205 

17.664 

6.870 

1519. 13 

5 1 . 208 

’  210 

17.678 

6.797 

1515. 36 

50.578 

215 

17.693 

6.710 

1511.61 

49.894 

220 

17.707 

6.611 

1507.99 

49.035 

225 

17.721 

6.498 

1 504 . 4 1 

48. 121 

230 

17.736 

6.373 

1 500 . 90 

47.125 

'77'cr 

17.748 

6.236 

1497.47 

46.  O  c*  7 

240 

17.761 

6.087 

1497.47 

44.870 

245 

17.774 

5.927 

1490. 89 

43.628 

250 

17.786 

5.  756 

1487.75 

42 . 308 

255 

17.798 

5.576 

1484.71 

40.929 

260 

17.810 

5.388 

1481.71 

39.498 

265 

17.821 

5.  192 

1478.95 

38 .012 

270 

17.832 

4 . 990 

1476.24 

36. 488 

275 

17.842 

4.783 

1473.65 

34 . 933 

280 

17.852 

4.573 

1471.18 

OO  i  9U«. 

60 


TABLE  3-1 

(CQNT ' Q) 

angle 

Re 

aht 

Tg-Tw 

• 

Q 

0 

0 

BTU 

it2 

°R 

BTU 

hr — -ft2— °R 

sec 

285 

17.861 

4.362 

1468.82 

3 1 . 786 

290 

17.870 

4.  152 

1466. 58 

30 . 226 

295 

17.879 

3.  944 

1464.46 

28.685 

300 

17.887 

3.  741 

1462.44 

27.198 

305 

17.895 

3.545 

1460.54 

25.737 

310 

1 7 . 902 

3.  356 

1458. 73 

24.346 

315 

17. 909 

3.  181 

1457. 03 

23. 057 

320 

17.915 

3.017 

1455.41 

21.851 

325 

17.921 

2.867 

1453.88 

20 . 750 

330 

17.927 

2.734 

1452.41 

19.774 

335 

17.933 

2.618 

1451.01 

18.923 

340 

17.938 

2 .  522 

1449.67 

18.217 

345 

17.944 

2.  445 

1448.36 

17.651 

350 

17.949 

2.  390 

1447. 09 

17. 244 

355 

17.954 

2.357 

1445.83 

16. 996 

TDC 

17.959 

2 . 346 

1444.58 

16. 906 
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(Q:BTu/sec) 
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CHAPTER  FOUR 


RADIANT  HEAT  TRANSFER  IN  THE  NEW  CYCLE  ENGINE 
EXPANSION  CYLINDER 

1 .  THEORETICAL  BACKGROUND 

When  a  temperature  difference  exists  between  two 
bodies,  part  of  the  thermal  energy  is  transferred  as 
electromagnetic  waves  of  radiation  at  varying  wavelengths,  A, 
characteristic  of  the  body  temperature  and  identity.  In  the 
el ectromagneti c  spectrum,  thermal  energy  is  propagated 
through  .01  >im  *=  A  ^  100  )un.  Bodies  above  absolute  zero 
temperature  continuously  emit  energy  which,  when  striking 
another  body,  may  be  partly  absorbed,  to,  reflected,  t,  or 
transmitted  ,  T,  where: 


to  +  ^  +  T  =  1  (39) 

Thermal  radiation  reflected  from  and  transmitted  through  a 
body  is  again  absorbed,  reflected  by,  or  transmitted  through 
any  additional  bodies  ultimately  establishing  an  equilibrium 
state.  An  ideal  body  that  absorbs  all  radiation  is  known  as 
a  black  body;  if  no  radiation  is  transmitted,  as  in  the  case 
for  most  industrial  liquids  and  materials,  the  body  is  known 
as  ooaque.  Combustion  chambers  covered  with  a  thin  film  of 
oil  contaminated  with  soot  particles  will  absorb  a 
significant  percentage  of  incident  radiation,  T  =  0. 

Thermal  energy  emitted  from  a  body  unit  surface  per 
unit  time  at  a  particular  wavelength,  the  monochromatic 


emissive  power, 


can  be  determined  as  a  function  of 


EX  ’ 

temperature  T(°K)  by  Planck's  Law: 


=  2ircc  2X  5Ccc„/ecbxT-i:]  1 


(40) 


erg /sec -cm"' 


where: 


=  speed  of  light  2.997902x10*^  cm/sec 

—  ^7 

=  Planck's  constant  6.623777x10  “  erg-sec 


cb  =  Boltzmann's  constant  1 . 38026x lO-1^  erg/°k 
X  =  wavelength  in  cm 

T  =  temperature  in  °K 

will  increase  with  increasing  X  to  a  maximum  value  then 
decreases,  the  distribution  shifting  toward  the  shorter 
wavelengths  at  higher  body  temperatures.  The  relationship 


between  T(°K)  and  X 


Displacement; 


ma;< 


is  determined  by  Wein's  Law  of 


Xmax  T<°K>  =  2898  cm-°K 

The  total  radiant  emissive  power,  E,  may  therefore  be 
calculated  at  a  given  body  temperature  by  integrating  the 
area  under  the  curve: 

E  =  f"  E.  dX  (41) 

Jo 

For  a  given  black  body,  substitution  of  equation  (40)  into 
equation  (41)  yields  the  Stephan  Boltzmann  Law: 

EBlack  Bod^  ”  ^o^  (°R)  (42) 

BTU/hr-f t2 

where:  <r  =  Stephan  Boltzmann  constant  1 . 713x  10-<?BTU/ht — f  t2°R4 


Real  bodies,  however,  are  not  absolutely  black.  Those 
bodies  of  a  continuous  radiating  electromagnetic  spectrum 
where  the  curve  is  proportional  to  that  o-f  an  absolute 

black  body  at  each  A  are  known  as  grey  bodies.  The  ratio  of 
emissive  power  o-f  a  grey  body  to  an  absolute  black  body  is 
known  as  the  coe-f -f  1  ci ent  of  emissivity,  £: 

£  grey  body  (T)  =  £  (43) 

E  black  body  <T) 


From  an  analysis  of  equilibrium  emittance  and  absorptance 
between  a  grey  body  and  its  surroundings  at  the  trivial  case 
of  each  equal  to  the  same  temperature,  it  can  be  demonstrated 
that  €  is  equal  to  Jo  (Kirchoff's  Law).  This  important 
phenomenon,  strictly  correct  only  within  the  defined 
limitations,  is  fortunately  ooserved  to  be  approx i matel y 
correct  in  many  real  conditions  where  temperature  differences 
are  not  severe.  Additionally,  many  real  industrial  materials 
radiate  near  grey  at  very  high  values  of  Jo;  instantaneous 
radiant  heat  transfer  rate,  GR  ,  between  two  near-grey  bodies 
1  and  2  of  different  temperatures  at  a  distance  apart  may 
therefore  be  calculated  from  equation  (44): 


Q 


Rl,2  =  G1  .Z^o^l  ,2AHT  <T14-T24>  ( 1  hr/,360°  sec) 

Btu/sec 

(44) 


■1,: 


1  +  1-1 


G1  e2 


F,  =  geometric  view 
hi 


factor  between  bodies  1  and  2 
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Tables  providing  -formulas  -for  £  and  F  calculation 


between  multiple  bodies  at  various  common  geometries  are 
published  in  standard  heat  transfer  texts.  For  gas  or 
incandescent  solid  particles  radiating  as  a  de-fined  unit  body 
within  an  encompassing  simple  geometry  combustion  cylinder,  F 
may  be  assumed  equal  to  unity. 

2.  NON— LUMINOUS  GAS  RADIATION 

2. 1  CHARACTERISTIC  HEAT  TRANSFER  EQUATIONS 

Unlike  solid  bodies  which  radiate  across  continuous 
wavelength  spectrums,  gases  emit  and  absorb  energy  at  direct 
wavelengths  only  in  line  spectra.  For  simple  monatomic  and 
diatomic  gases  such  as  O^N^,  and  H.-, ,  the  line  spectra  are 
negligible  and  may  be  considered  transparent  (T  =  1).  For 
complex  gases  such  as  generated  in  exhaust  products,  however, 
radiant  capacity  may  be  significant.  CO^  and  H-.D  vapor 
predominate  as  radiant  emitting  gases  in  combustion  products: 
the  significant  bands  of  radiation  are  C141: 


C02 _  _ h2o 


Band 

1 

2 . 36  —  3 . 02 

2.24  -  3.27 

>tm 

Band 

2 

4.01  -  4.80 

4.8  -  8.5 

>tm 

Band 

3 

12.5  -  16.5 

12  -  25 

>im 

In  radiating  gas,  both  emittance  and  absorptance  exist 
throughout  the  gas  volume.  Penetrating  radiant  energy  will  be 
reduced  by  absorption  proportional  to  the  number  of  molecules 
in  the  ray  path  as  a  direct  function  of  individual  gas 
partial  pressure,  temperature,  and  path  length,  X.  This 
intensity  attenuation  is  expressed  by  Beers  Law: 
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based  -fuels, 


the  resultant  total  gas  emissivity  may  be 


considered  equal  to: 


£ 


9 


(49) 


Again,  significant  data  base  has  allowed  development  of 

graphs  for  £  ,  £  ,  and  as  a  function  of  P,  Tq  and 

CCU  hUO  9 

jL 

lr,,  and  are  published  in  standard  heat  transfer  texts. 
Empirical  formulas  based  on  experimental  data  for  calculation 
of  emissive  power  of  CCU  and  l-UQ  vapor  useful  for  a  first 
approx i mati on  in  U.S.  units  are  1141: 

Ec0  =  1 . 7 125x 10~3 (Pqq  1R) l/3(Tg/100>3,5  <50) 

2  2 
BTU/sec-ft^- 

EH  0  =  1.6738xl0-5  PH  Q-8lR*6(Tg/100)3‘°  (51) 

2  n  2 

BTU/sec-f  f*" 

2.2  PREDICTION  OF  GA5  RADIATION  HEAT  TRANSFER 

References  C3,  6,  8,  9,  10,  14,  15,  163  state  that  gas  . 

radiation  is  insi gni f i cant  in  calculating  high  rc  Cl  engine 

heat  transfer  loss.  To  ascertain  if  separate  consideration 
of  gas  radiation  in  the  low  rc  New  Cycle  Engine  expansion 
cylinder  may  be  required,  it  is  sufficient  to  estimate  QR 

magnitude  at  maximum  1R  and  Tg  in  the  exhaust  stroke  and 

compare  to  Q  calculated  by  Woschni's  correlation.  For  a 
typical  diesel  fuel  of  C/H  mass  ratio  equal  to  6  operating  at 
lean  #,  distribution  of  combustion  product  gas  includes 
approximately  10'/.  CCU  and  107.  H20  water  vapor  C143. 


Calculation  o-f  QR  using  equation  (48)  and  (49)  at  Appendix  A 
parameter  values  at  exhaust  stroke  BDC  -follows: 

Q  Woschni  =  52.909  BTU/sec 
V  =  1.43333  ft3 

Aht  =  7.037  -ft2 

Tg  =  2439. 06° R 

P  =  1.05  atm:  F'H  q  =  .105  atm  PCq  =  .105  atm 

2  2 

1R  =  3.4(1.43333)  =  .6925  ft 

7.037 

Qco  «  ( 1 . 7125x 10~5) C (. 105) (.6925) 3 1 /3 , 2439 . 06 , 3‘ 5 ( 7 . 037 ) 

2  100 

=  3.604  BTU/sec 

□H  q  *  ( 1 . 6738x 10~5) (. 105) ‘ 8 (.6925) ' 6 ,2439. 06, 3' 0 ( 7. 037 ) 

2  100 

=  0.226  BTU/sec 

□R  «  3.830  BTU/sec 
gas 

In  the  above  calculations,  corrections  -for  ££g,  Ew,  and 
emissive  power  o-f  the  cylinder  surface  proportional  to  i^Tw4, 
all  which  would  reduce  the  total,  have  been  ignored.  The 
estimated  value  equal  to  3.83  BTU/sec  is  well  absorbed  into 
52.909  BTU/sec.  Noting  that  -for  all  other  crank  angles  in 
the  exhaust  stroke  both  Tg  an  1R  decrease,  plus  in  the 
combustion-expansion  stroke  luminous  -flame  radiation  exists, 
it  may  be  assumed  that  separate  consideration  o-f  gas 
radiation  in  the  New  Cycle  Engine  expansion  cylinder  heat 
transfer  calculation  is  not  required. 


LUMINOUS  FLAME  RADIATION 


3. 1  FUNDAMENTAL  CAUSE 

The  phenomenon  of  Flame  radiation  is  quite  complex, 
dependant  intimately  on  the  physical  and  chemical  properties 
oF  combustion,  pressure,  temperature,  and  a  myriad  oF  engine 
and  Fuel  parameter  variables.  Two  basic  Flame  types  are 

produced  with  respect  to  radiant  capacity:  non— luminous  Flame 
character  1 st i c  oF  SI  engines,  and  luminous  Flame 
character i sti c  oF  Cl  engine  diFFusion  Flames.  In  SI  engines, 
pre— mixed  Fuel  and  air  are  adi abati cal  1 y  heated  in 

compression  then  ignited  by  a  spark  resulting  in  combustion 
by  turbulent  propagation  oF  Flame  through  the  charge  From  the 
point  oF  ignition.  The  mixture  burned  at  the  beginning  oF 
Flame  travel  expands,  then  recompresses ,  reaching 
temperatures  hundreds  oF  degrees  higher  than  that  which  is 
burned  last.  This  recompressi on  oF  combustion  product  gases, 
particularly  C02,  yields  visible  yellow  radiation.  No 

intermediate  solid  particles  are  Formed,  however,  and 

apparent  emissivities  oF  combustion  product  gases  as  a 
Function  of  near  transparent  A  bands  is  small.  The  driving 

4 

potential  For  radiant  heat  transFer,  proportional  to  T  , 
diminishes  more  rapidly  then  Farced  convection  heat  transFer 
proportional  to  T;  the  contribution  oF  instantaneous  radiant 
heat  transFer,  Q^,  integrated  over  the  engine  cycle  is 
insigniFicant. 

In  Cl  engines,  the  Fuel  spray  injecting  into 
adi abatical ly  compressed  air  near  piston  TDC  is  ignited  by 
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spontaneous  explosion  at  a  mi cro- mi x ture  of  fuel  and 
entrained  air,  thereafter  burning  diffusively.  There  are 
many  sources  of  ignition  and  flame  tends  to  conform  to  the 
fuel  spray  shape  until  dispersed  by  gas  motion.  In  the 

chemical  combustion  of  hydrocarbon  based  fuels,  cracking  to 
simpler  structures  and  polymerization  occurs.  When 

insufficient  local  0.-*  exists  to  feed  the  rapidly  increasing 
reactivity  of  the  mixture,  thermal  decomposition  of  C-H  ana 
C-C  bonds  occurs  resulting  in  incandescent  carbon  soot 
formation  radiating  in  a  continuous  spectrum.  Alcock  et  al 
[111,  in  a  classic  high  speed  photographi cs  study  of  a  Meurer 
direct  injection,  medium  speed  Cl  engine,  clearly  observes  a 
luminous  flame  front  filling  the  chamber  within  20°  of 
ignition  start  <10°  BTDC) ,  glowing  white  hot  at  peak  flame 
temperatures  equal  to  4500° R ,  then  cooling  by  expansion  to 
dull  red.  Unlike  SI  engines  where  combustion  completes  early 
in  the  expansion  stroke,  the  Cl  combustion  process  extends 
late  into  the  expansion  stroke. 

This  difference  between  local  02  concentration  and 
cracking  rate  is  the  principle  cause  of  soot  formation  and 
luminous  flame  radiation.  The  magnitude  of  associated  heat 
transfer  cannot  be  ignored,  and  has  been  measured  by  Annand 
C33,  Sitkei  C143,  Dent  et  al  CIO},  Flynn  et  al  C153,  and 
Qguri  et  al  C83  as  typically  20'/.  to  307.  of  total  heat 
transfer  loss,  Qy  ,  in  non-swirl  enhanced  four  stroke,  direct 
injection,  high  rc  Cl  engines  under  varying  design  and 
operational  conditions.  Oguri  et  al  C81  demonstrate  this  in 


F i gure  4- 1 , 


showing  QR/Qy  between  two  such  different  size  Cl 
engines  operating  at  varying  RF'M ,  load,  and  cylinder  wall 
soot  contamination:  Engine  A  is  a  5.9  inch  bore  natural 

aspirating  engine;  Engine  B  is  a  12.2  inch  bore  supercharged 
engine.  In  very  large  diesel  engines  the  contribution  to  Qy 
of  radiant  heat  transfer  is  even  greater;  J.  Butler  of 
Dexford  Co.,  United  Kingdom,  estimated  the  radiant  component 
o-f  Qy  as  75V.  -for  a  40  inch  bore  diesel  C83.  In  any  case, 
measured  peak  -flame  emi  ssi  vi  ti  es ,  £y,  approaches  unity 

(1.0),  and  while  reaction  zone  temperature  is  di-f-ficult  to 
determine,  a  consensus  opinion  shared  is  that  peak  radiant 
temperature  is  much  higher  than  Tg.  Chapman  et  al  state 
C17,  pg  133: 

"the  apparent  radiant  temperature  is  much  closer 
to  the  flame  temperature  than  the  average  bulk 
gas  temperature. " 

Peak  flame  temperature,  Ty  ,  reached  in  high  r_  Cl  engines  is 
typically  4200  ±  200° R. 

3. 2  KHAN  ET  AL  COMBUSTION  SOOT  PRODUCTION  MODEL 

While  a  valid  quantitative  calculation  of  soot 
concentration  in  New  Cycle  Engine  combustion  is  not  possible, 
a  brief  study  of  a  soot  production  model  for  a  simple  design, 
direct  injection,  quiescent  Cl  engine  will  illustrate  the 
critical  parameters  on  which  soot  formation  is  based  and 
establish  a  foundation  for  empirical  estimation  of  peak  New 
Cycle  Engine  Q^.  Khan  et  al  C183  have  developed 


semi -empirical  equations  for  soot  formation  in  a  fuel  jet 
mixing  model  applied  successfully  to  a  9.8  cm  bore  engine  at 


cl/  a 


FIGURE  4-1 s  Qr/ Qt  f°r  two  4  stroke,  direct  injection  Cl 

engines  under  varying  operating  conditions  as 
observed  by  Oguri  et  al  C8D 
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2000  RF'M.  This  simple  model  ignores  complex  stratification 
and  consideration  of  fuel  droplet  evaporation;  all  unburnt 
fuel  is  assumed  to  be  in  the  vapor  state.  The  jet  mixing  model 
injecting  into  a  cross-stream  air  flow  is  shown  in  Figure 
4-2.  The  air  entrainment  zone  (e)  is  c'omposed  of  a  fuel  rich 
zone  (c)  where  soot  formation  occurs,  a  products  zone  (d) 
where  intense  reaction  occurs,  and  an  outer  zone  composed 
primarily  of  air.  Outside  the  entrainment  zone  is  the  air 
mass  zone  (a) .  For  the  purpose  of  soot  formation 

cal cul ati ons ,  a  local  equivalence  ratio  of  unburnt  fuel  in 
the  fuel  rich  zone,  $  ,  is  defined.  A  series  of  equations 
are  proposed  by  Khan  et  al  for  #u  calculation;  the 
significant  points  of  which  is  that  $u  is  inversely 

proportional  to  air  entrainment  which  is  in  turn  a 
progressive  function  of  air  swirl  ratio  and  engine  speed.  §u 
is  not,  however,  directly  coupled  to  §  charge.  The  overall 
mechanism  of  soot  formation  under  high  temperature  and  short 
reaction  time  conditions  characteristic  of  diffusion  flame 
combustion  may  be  characterized  by  an  Arrhenius  type 
equation: 


d.3 

dt 


(52) 


where:d_S  =  soot  formation  rate 
dt 

Vu  =  volume  of  soot  formation  zone  (c> 

V  =  cylinder  volume 

F‘u  =  partial  pressure  of  unburnt  fuel  in  soot 
formation  zone  <c) 

Tu  =  local  temperature  of  soot  formation  zone  (c) 
cg  =  soot  formation  rate  coefficient 
Es  =  soot  combustion  activation  energy 
RQ  =  universal  gas  constant 
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FIGUR 


S  4-2:  Soot  production  fuel  jet 
Khan  et  al  Cl  3D 


mixing  model  of 


After  formation,  soot  particles  coagulate  and/or  combust  from 
a  typical  initial  diameter  of  250  A  where  combustion  of  soot 
coagulate  is  again  controlled  by  an  Arrhenius  type  equation 
proportional  to  the  partial  pressure  of  CU. 

The  values  of  constants  and  coefficients  are,  of 
course,  a  function  of  individual  engine  design  and  fuel 
structure  selected  to  best  "curve  fit"  experimental  data. 
Valuable  qualitative  conclusions,  however,  stand  out  from 
the  work  by  Khan  et  al : 

1.  The  rate  of  soot  formation  depends  on  three 

principle  parameters:  local  F’  (partial), 

T,  and  f  in  the  fuel  rich  soot  formation 
cone.  The  effect  of  varying  overall  $ 
charge  is  moderated  as'  §u  is  always  near 
unity.  Swirl  ratio,  however,  has  a 
significant  effect.  A  comparison  of  soot 
concentrations  measured  by  Khan  et  al , 

Flynn  et  al ,  and  Dent  et  al  is  shown  in 
Figure  4-3  which  graphically  illustrates 
soot  concentration,  thus  luminous  flame 
radiation,  suppression  in  high  swirl  Cl 
engines. 

2.  Approximately  607.,  as  measured  by  Khan  et  al 
C183,  of  soot  initially  formed  is  burnt  as 
it  moves  from  the  cooler  fuel  rich  zone  (c) 
into  the  products  zone  (d)  where  soot 
oxidation  incandescent  radiant  emission  is 
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a  grey  body,  optically  thick.  The  exhausted 
soot  is  that  portion  which  escapes  the  high 
temperature  products  reaction  zone  (d)  in 
the  process  of  turbulent  diffusion  from  the 
fuel  rich  to  the  lean  zone. 

3.  High  speed  photography  indicates  that  the 
majority  of  soot  formation  occurs  at  the 
beginning  of  the  main  combustion  phase, 
within  an  extremely  short  period  Csee 
Figure  4-31.  The  total  amount  of  soot 
remains  constant  in  the  later  combustion 
stages. 

A  clear  conclusion  of  this  model  is  that  both  soot 
formation  rate  and  soot  combustion  rate  are  progressive 
functions  of  fuel  and  02  partial  pressures.  The  New  Cycle 
Engine,  operating  at  signif icantl y  reduced  expansion  cylinder 
pressures,  should  therefore  experience  a  marked  reduction  in 
luminous  flame  radiant  heat  transfer.  Due  to  lack  of  (built) 
engine  test  data  however,  the  degree  to  which  QR  decreases 
relative  to  magnitudes  observed  in  high  rc  Cl  engines  must  be 
estimated  from  indirect  empirical  rel at i onshi ps .  The  quality 
of  this  estimation,  based  principly  on  reduced  pressure 
consideration  in  low  swirl  conditions,  must  be  tempered  with 
the  knowledge  that  many  other  operating  parameters  have  a 
secondary  effect  on  flame  emissivity  and  temperature. 


o 


o  FLYNN  ET  AL 


•  DENT  ET  AL 
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ylinder  soot  conentration , 
gle  as  observed  by  various 


3.3  SECONDARY  FACTORS  AFFECTING  CQNBUST ION  SOOT 

CONCENTRATION 

Extensive  tests  performed  by  Flynn  et  el  on  a  four 
stroke,  direct  injection  Cl  engine,  with  piston  cup 
combustion  chamber  enlarged  to  suppress  swirl  ratio  and 
enhance  soot  production,  illuminate  many  of  these  secondary 
effects;  the  important  ones  of  which  are  listed  below  C161: 

1.  QR  increased  with  inlet  manifold  pressure 
increase  at  nearly  proportional  rate; 
much  less  ignition  delay  and  fuel /air 
pre-mixing  occurs  causing  increased 
tendency  for  soot  formation  in  diffusion 
f 1 ame. 

2.  As  injection  time  was  advanced  an 
increase  in  Tg  and  an  increase  in  peak 
apparent  radiant  temperature  occurred; 
due  to  fourth  power  effect  a  significant 
increase  in  peak  radiant  emissive  power 
resul ts. 

3.  Decreasing  cetane  number  with  fixed 
reference  fuel  families  encouraged 
ignition  delay,  increased  fuel/air  pre¬ 
mixing,  and  reduced  soot  formation;  at 
supercharged  operation  this  effect  is 
effectively  suppressed  due  to  predominant 
effect  of  manifold  pressure  on  ignition 
delay. 
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4.  Increasing  C/H  mass  ratio  fuels  at  fixed 
cetane  numbers  encouraged  UR  due  to 
simple  availability  of  carbon 

concentration;  the  rate  is  not 

proportional  to  C/H  mass  ratio  and 
apparent  only  at  high  ratio  aromatic 
fuel s . 

Flynn  et  al  additionally  note  that  the  apparent  flame  optical 
thickness,  k^l^,,  diminishes  approximately  in  inverse 
proportion  to  increasing  cylinder  volume.  This  conforms  with 
the  Khan  et  al  model  and  Figure  4-3,  indicating  that  the 
number  of  incandescent  radiating  particles  in  later 
combustion  stages  is  essentially  constant. 

3.4  CHARACTERISTIC  HEAT  TRANSFER  EQUATIONS 

Emission  as  a  function  of  wavelength  for  clouds  o+ 
small  particles  obeys  a  variation  described  by  equation  (53): 

Ea  '=  1  -  exp  <-kxlR/,\*95)  (53) 

The  spectrum  is  near  grey,  and  for  simplicity  of  calculation 
can  be  converted  to  equivalent  grey  body  emissivity: 

eAqAdx 

-  — i -  04) 

f*n 

j 

J*o 

where:  ?«0  -  ,\n  is  range  of  significant  A  band 

The  effect  of  combustion  cylinder  wall  reflectivity,  ?w,  as 

demonstrated  by  Qguri  et  al  in  Figure  4-1,  may  be  corrected 
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■for  by  determination  of  an  apparent  flame  emissivity 


accoraing  to  equation  (55): 


<1-Vef  -Mr 

£  =  -  =  1  -  e  (55) 

‘-Vf 

where:  k^  =  apparent  flame  absorption  coefficient 

per  unit  length 

When  flame  envelopes  the  entire  cylinder  volume,  observed  by 
Alcock  et  al  [11]  to  occur  approximately  10  to  20°  after 
ignition  start,  1R  may  be  calculated  by  equation  (47)  as  a 
function  of  instantaneous  cylinder  0  and  Sitkei  et  al 

C153  state  that  surfaces  of  pre-swirl  chambers,  piston  crown 
and  cylinder  heads  have  ?  =  0.1  (approx.);  cylinder  sleeve 

?w  =  0.2  (approx.).  It  should  be  noted  that  estimates  for 
vary  greatly,  in  part  due  to  different  modeling  of  cylinder 
surfaces  as  specular  reflective  to  semi-rough  diffusive 
reflective  characteristic  of  oil  soaked,  soot-contaminated 
surfaces.  QR  may  then  be  determined  by  equation  (56): 

Qr  =  £a^0^HTCTf  4-Tw4;5  (1  hr/3600  sec)  (56) 

^BTU/sec 


3.5  SITKEI  ET  AL  METHOD  OF  APPARENT  FLAME  EMISSIVITY,  g ^ , 
CALCULATION 

A  useful  empirical  relationship  used  by  Annand  C33  and 
Dent  et  al  C103  based  on  known  soot  concentration  is  stated 
in  equation  (57)  in  U.S.  units: 

kf  =  5898  &  ft"1  (57) 

where:  A  =  soot  concentrations  in  lbm/ft"* 
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For  estimation  oF  Sa  without  knowledge  o+  soot 
concentrati  on ,  Sitkei  et  al  propose  the  -following  definition 
o-f  Ea  as  a  -function  o-f  instantaneous  P,  #,  lp  and  crank  angle 
past  ignition  start  C153: 

- .  3048k  F'l  p 

£_  =  £_  Cl  -  e  P  3  <  58 ) 

a  o 

where:  SQ  =  emissivity  o-f  in-finitely  thick  radiating 
Flame;  ss  0.8 

kp  =  apparent  Flame  absorption  coeFf lcient  per  unit 
length-unit  pressure  (atm) 

Values  oF  kp  are  determined  From  Figure  4-4,  at  applicable 
crank  angle,  and  excess  air  Factor  equivalent  to  |-1.  Sitkei 
et  al  based  Figure  4-4  on  measurements  conducted  by  a  photo¬ 
electric  pyrometer  unit  through  a  quarts  view  window  insert 
to  a  swirl  pre-combustion  chamber  oF  a  Four  stroke  Steyr 
diesel  at  rc  equal  to  18.5.  Recognizing  that  £a  is  dependant 
on  many  additional  variations  in  engine  operating  parameters, 
Sitkei  states  that  to  a  First  approximation  C14,  pg.  563: 

"in  the  absence  oF  other  data,  it  is  Felt  that  k 
given  here  is  universally  adaptable  to  all 
engines  of  the  diesel  type....  with  the  known  P-V 
diagram  a*  the  engine  under  test  the  variation  oF 
kp  can  be  Found" 

Inherent  in  this  assumption  is  an  assumed  %  oF  107.  to  20/1, 
an  ignition  start  near  TDC,  a  rather  arbitrary  value  oF  £Q  = 
0.8  based  on  observed  steady  state  diesel  Fuel  Flame 
emissivities  in  combustors,  and  the  belieF  that  £a  observed 
in  a  2.6  cm  pre-combustion  chamber  may  be  extrapolated  to 
main  cylinder  emissivities. 


It  should  be  noted  that  the  values  o-f  determined  by 
the  above  method  are  generally  less  than  observed  by  a  number 
o-f  other  investigators  in  direct  injection  Cl  engine  testing. 
Dent  et  al  state  CIO,  p.  17603: 

"The  essential  results  o-f  Sitkei  and  Ramanaiah 
are  in  agreement  with  those  o-f  Flynn  et  al  in 
terms  of  the  trend  of  flame  emissivity  and 
optical  thickness  with  crank  angle,  load,  and 
overall  air  fuel  ratios.  However,  the  magnitude 
of  the  flame  emissivity  and  optical  thickness  are 
approx i matel y  half." 

For  first  approximation  of  predicted  peak  £a  in  the  New  Cycle 
Engine  expansion  cylinder,  it  is  therefore  believed  that  the 
£a  obtained  by  the  Sitkei  et  al  method  should  be  multiplied 
by  a  factor  of  2.  It  is  also  noted  that  Sitkei  et  al  agree 
with  typical  QR/Qy  range  of  207.  to  307.  even  at  reduced  £a 
values  due  to  assumption  of  Belinskiy's  significantly  higher 
T+  profile  shown  in  Figure"  4-5,  vice  direct  optical 
measurement . 

3.6  LOW  COMBUSTION  PRESSURE  POWER  UNIT  METHOD  OF 
APPARENT  FLAME  EMISSIVITY,  £  ^ ,  CALCULATION 
For  first  approximation  purposes,  an  alternate  method 
of  predicting  £a  is  to  examine  the  characteristic  emissivity 
equations  used  in  power  units  operating  at  the  low  peak, 
combustion  pressures  of  the  New  Cycle  Engine.  Different 
properties,  continuous  flaming  and  non-reciprocating 
mechanics  of  gas  turbines  do  not  allow  broad  application  of 
flame  emission  data  to  internal  combustion  engines.  However, 
actual  combustion  kinematics  and  flame  luminosity  observed  in 
power  units  at  pressure  conditions  more  reflective  of  the  New 
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Cycle  Engine  should  not  be  ignored.  It  is  therefore 
considered  that  a  -fair  prediction  of  New  Cycle  Engine  peak  Ea 
only  can  be  ascertained  by  use  of  gas  turbine  power  unit 
emissivity  equations  using  diesel  ■fuel  properties. 

For  gas  turbine  application,  £a  can  be  determined  by 
equation  <59)  C19D: 


£ 


a 


where:  z 


1  -  e  z  (59) 

-38636.55  PC  (fuel/air  mass  r  at  i  o)  1  ^,3  *  JTg-*  * 


The  flame  luminosity,  Z,  can  be  determined  as  a  function  of 
fuel  C/H  mass  ratio  from  Figure  4-6.  For  commercial  liquid 
fuels  of  conventional  character  ranging  from  gasoline  to  fuel 
oil,  C/H  mass  ratio  can  be  estimated  using  equation  (60): 


C  mass  ratio  =  . 16015p  -  1.66  (60) 

H 

where:  p  fuel  density  in  lbrn/ft*'  is  referenced  to 

518°R  (288°K) 

3. 7  BOLTZMANN  NUMBER  FUNCTION  METHOD  OF  FLAME 
Temperature,  Tf ,  CALCULATION 

A  consensus  opinion  of  i nvesti gators  support  that  for 
first  approximation  purposes,  radiant  temperature  is  equal  to 
Tjj  various  T^  profiles  published  from  high  rc  Cl  engine 
testing  are  shown  in  Figure  4-5.  An  empirical  correlation 
allowing  calculation  of  T^  is  the  Boltzmann  number  function, 
Tg/T^  =  f  (Bo*) ,  established  in  Figure  4-7.  The  Bo*  is 
defined  in  accordance  with  equation  (61)  E143: 
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Bo*  = 


^e^p  -fuel 


roc'-fMHT  1  f 


where: 


B  =  -fuel  consumption  rate 


C  =  mean  total  specific  heat  of  combustion 

products  o-f  1  lbm  o-f  -fuel 

For  ease  of  calculation,  B  and  C  -fuel  may  be  calculated  -for 

e  p 

the  New  Cycle  Engine  using  the  -following  identities: 


lbm/hr  -fuel 


(SFC)  (.Aheat)  <mcharge)  (NRpM)  <60  min /hr) 
2544  EtTU/hr-HP 


C„  fuel  at  (Air/Fuel  mass  ratio)  (C„  combustion  products)  (63) 
p  P 

Substitution  of  the  above  identities  into  equation  (61) 
yields: 


^  C(SFC)  (Aheat)  (mcharge)  (Np,pM>  (Air/fuel  mass  ratio)Cp_ 


Bo  = 


where:  C  =  1.3768x10 


£f AHTTf " 


_  +4 

Calculations  of  at  a  known  Tg  and  from  the  Bo  function 
requires  an  iterative  process  matching  an  assumed  T+  both 
directly  from  the  generated  Tg/T^  fraction  and  its  coupled 
Bo  from  Figure  4-7.  An  increasing  E^  will  decrease  the  Bo  , 
thus  decrease  Tg/T+,  resulting  in  an  increased  calculated  . 
It  is  noted,  however,  that  increases  in  E^  and  are  not 
proporti onal :  over  a  wide  range  of  Bo*  (  >  .  6  ),  Tg/T+  slope 

M 

is  only  slightly  positive.  At  Bo  greater  than  1.2, 
characteristic  of  low  luminous  flames,  the  value  of  is 
essentially  independent  of  Bo*  (thus  £+)  and  equal  to  Tg/0.8. 
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3.8  PREDICTION  OF  PEAK  COMBUSTION  £ ,  and  T, 

For  calculation  of  peak  £a  and  T^  in  the  New  Cycle 
Engine  expansion  cylinder  at  -full  load,  the  following 
assumptions  are  made: 

1.  The  expansion  cylinder  ref  1 ecti vi ty ,  T  =  0.1 

2.  The  flame  growth  encompasses  the  entire 

cylinder  volume  at  20°  maximum  crank  angle  from 
ignition  start  Ccrank  angle  =  65° 3;  flame 

growth  rate  is  constant. 

3.  The  flame  luminosity,  £,  calculation  is  based 
on  medium  speed  diesel  fuel  C/H  mass  ratio  =  6.3 

4.  The  fuel /air  mass  ratio  =  .05  (1/20)  for  $  =  . 75 

5.  C_  at  fuel  C/H  mass  ratio  =  6.3  may  be 

Pep 

approximated  by  equation  (32)  as  a  function  of  Tg 
for  CnHVn  fuels  (C/H  mass  ratio  =  6.0). 

The  Sitkei  et  al  method  of  £a  calculation  was  developed  on  a 
test  engine  exhibiting  typical  high  rc  combustion  pressure 
profiles,  and  TDC  ignition  start,  both  conditions  not  met  in 
New  Cycle  Engine  operation.  The  peak  kp  from  Figure  4-4  at 
the  correct  excess  air  factor  va) ue  will  therefore  be  applied 
at  the  crank  angle  of  maximum  combustion  pressure  where  flame 
growth  reaches  100*/.  cylinder  volume. 

Sitkei  Method;  E_  peak 

*£a  peak  =  *51 

excess  air  factor  =  (.75)-1  =  1 . 33 

peak  kp  at  1.33  excess  air  factor  =  0. 63 

(see  Figure  4-4) 


Using  equations  (47)  and  (58),  plus  applying  a  Ea  correction 
■factor  of  2: 


•  *  V  flHT  1„  kp  .3048kpPlR  Ja  2£, 

at*  ft3  ft2  ft  (ft-at«fl 

45*  3.76064  .489331  3.944  .4210  .63  .304  ,210  .420 

peak  =  .420 
at  9  =  65° 


Low  Combustion  Pressure  Power  Unit  Method,  £  peak 

■flame  luminosity,  <L  =  4 . 0 
(see  Figure  4-6) 


Using  Equations  (47)  and  (59): 


9 

P 

ati 

Tq 

•R 

V 

ft3 

A 

ft2 

flaie  growth 

1 

V* 

ft 

l 

£a 

30 

4.24684 

2259.6 

.309053 

3.358 

25 

.0782 

.382 

.318 

35 

4.09172 

2371.21 

.366348 

3.345 

50 

.1757 

.453 

.366 

60 

3.92401 

2864.13 

,426316 

3.741 

75 

.2906 

.477 

.379 

65 

3.76064 

3136.02 

.488331 

3.944 

100 

.4210 

.480 

.381 

70 

3.26107 

3064.15 

.351767 

4.152 

too 

.4518 

.447 

.360 

peak  £a  =  .381 
at  0  =  65® 

A  repeated  calculation  -for  both  methods  at  peak  combustion 
pressure  greater  than  50  atm,  typical  -for  direct  injection, 
high  rc  Cl  engines,  results  in  peak  £a  x  1.0,  which  closely 
agrees  with  the  expected  value  (:*0.9)  . 

Boltzmann  Number  Function  Method,  peak 

±f 

The  known  values  required  for  calculation  of  the  Bo  at  65® 
crank  angle,  peak  £a  of  (.420  +  .381)/2  =  .40  are: 
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1 


4. 


Tg  =  3136.0  °R 

SFC  =  .284298  lbm  fuel/HP-hr 
Aheat  =  500  BTU/lbm  charge 
mcharge  =  -^24537  lbm 
NRPM  =  813-3 

C_  =  . 247+. 243x 10~4 (3136.0° R)  =  .3232  BTU 

PcP  lbm— ®"ft 


7. 

8. 


air  /-fuel  mass  ratio  =  20 


Aht  =  3.944  ftJ 


Substituting  the  above  values  into  equation  (64)  yields: 


Bo*  = 


(1.3768xl07) (.28429) (500) (.024537) (818.3) (20) (.3232) 
<3.944)£fTf3 


6. 435x10 10 
£fTf3 

From  equation  (55)  ,  an  estimated  peak  Ea  of  0.40  at  an 

assumed  =  0.1  yields  a  peak  equal  to  0.426.  The 

predicted  peak  may  therefore  be  calculated  -from  iteration 

£ 

using  Figure  4-7  at  Bo  equal  to: 


Bo*  = 


6.435xl010  =  l.SUxlO11 


.426T+-> 


V 


Assuming  a  value  o-f  =  Tg/0.8  =  3920° R,  the  calculated  Bow 

H 

equal  to  2.51.  This  Bo  is  greater  than  1.2,  therefore  the 
estimated  peak  T^  equals  3920°R.  Rounding  of-f  for  first 
approximation: 

peak  T^  =  4000° R 
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To  check  the  accuracy  o-f  this  procedure,  a  peak  £ ^ 
value  o-f  1.0  is  assumed  commensurate  with  observed  peak  in 
direct  injection,  high  rc  Cl  engines; 

6.435xl010  =  6.435xl010  . 

Bo#  =  - _  - - - 

(1.0)Tf3  T^3 

From  Figure  4-5,  the  observed  peak  is  typically  equal  to 
4200°R;  Tg/T^  now  equals  0.75.  Substituting  4200* R  into  Bo^ 
=  6.  435x  lO^/T^ yields  a  value  o-f  0.87  which,  entering 

Figure  4-7,  extrapolates  to  Tg/T+  =  0.76.  The  iteration 

essentially  matches,  indicating  that  a  peak  £,  and  peak  T, 
equal  to  1.0  and  4200° R  respectively,  are  in  -fact  reasonable 
values  -for  direct  injection,  high  rc  Cl  engines. 
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4. 


CALCULATION  OF  ORDER  OF  MAGNITUDE  REDUCTION 


IN  COMBUSTION  RADIANT  HEAT  FLUX 

Waschni  did  not  separately  measure  radiant  heat  flux  in 
his  engine  testing;  for  ease  of  calculation  he  simply 
absorbed  the  radiant  component  into  the  forced  convection 
heat  transfer  model  by  use  of  time  averaged,  "curve  fitting" 
multiplicity  coefficients.  Woschni  in  fact  believed  that 
C9 ,  pp.  3074  -  30751 s 


"the  total  radiation  amounts  to  no  more  than  207. 
of  the  total  combustion  term. .. In  steady  flames 
there  are  usually  large  combustion  chambers  with 
great  volumes  of  flames,  and  because  of  the 
small  velocities  there  are  small  convective  heat 
transfer  rates  only,  so  that  in  these  cases 
radiation  is  responsible  for  an  appreciable 
amount  of  the  total  heat.  Because  of  the  high 
velocity  and  density  of  gas,  the  convective  heat 
transfer  coefficients  in  internal  combustion 
engines  are  about  ten  times  higher,  so  that 
despite  the  higher  temperatures  the  radiation 
plays  an  appreciably  smaller  part. " 

While  the  thrust  of  the  statement  is  correct,  Woschni 
unfortunately  underestimated  QR ,  basing  cursory  radiant  heat 
transfer  calculations  on  the  assumption  of  £a  =  0.6  from 
steady  state  flame  observations  and  employing  for  the  unknown 
radiant  temperature,  Tg,  vice  T^. 

For  the  purpose  of  this  calculation,  establishing 
Woshni  peak  £a  and  values  of  1.0  and  4200° R  better 
conforms  to  published  data  on  Cl  engine  testing  of  similar 
design,  direct  injection,  high  rc  character i sti cs.  It  is 
noted  that  actual  observed  peak  Ea  values  are  typically  equal 
to  0.9  vice  1.0.  To  maintain  continuity  with  £fi  prediction 
methods  used  in  section  3.B,  however,  a  value  of  unity  is 
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assumed  as  predicted  for  high  peak  pressure  conditions. 


Without  knowledge  of  the  radiant  component  of  qw  and  QNqE  as 

a  function  of  crank  angle,  it  is  not  practical  to  attempt 

calculation  o-f  a  mean  order  o-f  magnitude  reduction  ratio, 

qW/,qNCE’  over  the  effective  90°  crank  angle  domain  where 

luminous  radiant  heat  transfer  can  be  significant  (£  a  0.1) 

a 

in  high  rc,  medium  speed  Cl  engines.  Fortunately,  in  this 
case  it  is  not  required.  Making  the  reasonable  assumption  of 
equivalent  £a,  decay  rates,  flame  emissivity  in  the  Mew 

Cycle  Engine  expansion  cylinder  should  certainly  decay  below 
0.1  at  advancing  crank  angles  from  ignition  start  prior  to 
+90°,  Thus,  peak  q^qNCE  should  be  a  minimum ,  conservative 


estimation  of  the  mean.  Therefore,  from  equation  (b6)  : 


qw 

—  «  —  peak: 

qNCE  qNCE 


« 


£a (Tf 4“^w4) NCE 


peak 

radiation  T  =  900° R 

-  w 

Substituting  in  the  respective  Ea  and  T^  values  previously 
determi ned : 


qw  1 . 0  <  42004— 9004 ) 

minimum  -  —  -  =  3.047  approx . 

qNCE  °*  4<40004-9004) 


AD-A144  831 


UNCLASSIFIED 


PREDICTED  HEAT  TRANSFER  LOSS  IN  THE  EXPANSION  CYLINDER  0/p. 
OF  A  TWO  CYLINDER. . (U)  MASSACHUSETTS  INST  OF  TECH  ^ 

CAMBRIDGE  DEPT  OF  OCEAN  ENGINEERS..  P  C  JORGENSEN 
MAY  84  N66314-70-A-0073  F/G  21/5  NL 


■ 

- - 

END  ! 

a -84 

DUO 

CHAPTER  FIVE 


CALCULATION  OF  CORRECTED  HEAT  TRANSFER  LDSS/HEAT 
RELEASE  PER  CYCLE,  Tt* 

Total  heat  transfer  loss,  Q-j-,  for  a  natural 
aspirating,  four  stroke  Cl  engine  includes  heat  transfer 
during  the  induction  and  compression  strokes  which  do  not 
exist  in  the  two  stroke  New  Cycle  Engine  expansion  cylinder 
cycle-  In  the  induction  stroke,  the  heat  transfer  driving 
potential,  Tg-Tw,  is  negative;  induction  charge  heating  is 
boosted  by  residual  gas  mixing.  In  the  compression  stroke, 
adiabatic  compression  may  boost  Tg  to  greater  than  1500°R. 
From  Figure  2-1,  it  is  observed  that  a  significant  positive 
Tg-Tw  occurs  only  in  the  last  30°  of  the  compression 
stroke.  The  order  of  magnitude  of  this  positive  driving 
potential  compared  to  that  observed  in  -the  combustion- 
expansion  and  exhaust  strokes  is  still  quite  small.  It  may 
therefore  be  assumed  that  the  induction  and  compression 
stroke  contribution  to  Q-j-  for  a  natural  aspirating,  four 
stroke,  high  rc  Cl  engine  is  not  significant.  A  direct 
coupling  of  observed  QR/Qy  ratios  may  then  be  applied  to  the 
New  Cycle  Engine  expansion  cylinder. 

From  Chapters  2,  3,  and  4,  the  following  principal 

values  have  been  determined: 

1.  fly  uncorrected  =  19.7617.;  Woschni  correlations 

2.  Qp/Qy  Woschni  engine  s  V;  .2  &  V  &  .3 


T 


3.  Forced  convection  qw/qNCE  =  6.087 

4.  Radiation  A^^NCE  =  minimum 

Uj  must  be  corrected  -for  the  -failure  o-f  Woschm's 
correlation,  based  solely  on  a  theoretical  -forced  convective 
heat  transfer  model,  to  properly  scale  down  luminous 
radiation  heat  transfer  to  low  rQ  operating  conditions. 

Correction  to  Tly  =  Tlp^  +  Tip 
At  high  r operation; 


qr 

QR/(Aheat) (mcharge) 

*R 

=  ¥ 

qt 

(Gr  +  QFC)/(^heat) (mcharge) 

ftR  +  aFC 

tlR  =  ¥TlR 

+  WFC 

TIr(  l-¥> 

"  VtlFC 

-  -SL 

rfc 

(65) 

Substituting  into  equation  (65)  the  known  order  of  magnitude 

reduction  ratios  for  low  rc  operating  conditions  and 

establishing  equality: 

At  low  r  operation: 

3.047  ftR*  =  _¥_  6.087  flFC* 
l-¥ 

ftR*  =  2.00  ¥  nFC*  (66) 

maximum  l-¥ 

As  Woschni's  correlation  properly  scales  down  forced 
convection  heat  transfer  only,  tlRC*  =  (1  -  ¥)fty.  Therefore  a 
conservative,  first  approximation  prediction  of  corrected 
total  heat  transfer  loss/heat  release  per  cycle,  tlT*, 


for  .2  6  V  6  .3  is 


1 


Farced  convection  heat  transfer  loss,  ftFC* 
RFC*  ~  <1-V)TIt  =  <1-V)  19.7617. 

13.837V.  £  tlFC*  £  15.809'/. 


Radiant  heat  transfer 


ftp.*  =  2. 00  <  V  ) ftFC* 
1  -  V 


loss,  UR* 

2.  QQ(  V  )  <  1-W)  tlT 
1-V 


=  2.00  V  tlT 

=  2.00  V  19.761V. 

7.904V.  £  ftR*  £  11.857V. 

max i mum 

Total  heat  transfer  loss,  ftT* 

ftT*  =  ftR*  +  ftFC*  =  2. OOVftj  +  (l-W)ftT 

=  <2. 00W+1— V) ftT 

=  <l+W)ftT 

=  (1  +  V)  19.761% 

23.713V.  6  ftT*  ~  25.689V. 

max i mum 
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CHAPTER  SIX 


CONCLUSIONS  AND  RECQMMENDAT I  DIMS  FOR  FUTURE  STUDY 

A  semi  -empi  ri  cal  correlation  -for  mean  instantaneous 
heat  transfer  coefficient,  developed  by  Woschni  on 
conventional  -four  stroke,  direct  injection,  quiescent  Cl 
engines,  was  modified  to  predict  heat  transfer  loss  in  a 
computer-generated  preliminary  design  engine  cycle.  The 
engine  divides  the  -four  stroke  -functions  between  a 
compression  cylinder  and  an  expansion  cylinder  interconnected 
by  a  static  regenerative  heat  exchanger.  Due  to  regenerator 
boost  heating  o-f  the  air  charge,  required  compression  ratios 
in  the  New  Cycle  Engine  expansion  cylinder  will  be 
si  gni  f  i  cantl  y  reduced,  and  increased  thermal  e-f-ficiency 
achieved.  To  verity  applicability  o-f  Woschni 's  correlation 
to  expansion  cylinder  heat  transfer ,  points  in  the  cycle  were 
examined  comparing  Woschni 's  correlation  to  -fully  developed, 
turbulent  -flow,  -forced  convection  relationships  -for  -flat 
plate  and  pipe  -flow.  The  mean  instantaneous  gas  velocity  was 
estimated  -from  the  construction  o-f  the  charge  impulse 
transient  existing  into  the  expansion  cylinder  -from  the 
regenerator  port.  An  order  o-f  magnitude  reduction  in  -forced 
convection  and  luminous  radiant  heat  -flux  was  calculated  to 
investigate  the  ability  o-f  Woschni  s  correlation  to  scale 
down  correctly  total  heat  transfer  to  low  compression  ratio 


operation. 


It  was  -found  that  the  impulse  velocity  transient  was 
sufficiently  suppressed  to  approximate  quiescent  conditions. 
It  was  found  that  laws  of  similarity  incorporated  into 
Woschni's  correlation  appear  to  allow  uncorrected  application 
at  low  compression  ratio  operation  for  forced  convection  heat 
transfer  only.  The  order  of  magnitude  reduction  in 
combusti on-expansi on  stroke  forced  convection  heat  flux  was 
found  to  be  as  much  as  twice  the  order  of  magnitude  reduction 
in  radiant  heat  flux.  A  correction  was  proposed  for  the 
application  of  Woschni's  correlation  to  achieve  a  conservative 
prediction  of  heat  transfer  loss  in  the  New  Cycle  engine 
expansion  cylinder  cycle: 

Qt  =  1.25  Qt  (67) 

corrected  uncorrected 

where:  Q-j-  =  T  Q  d$ 

uncorrected  J 

Q  =  he  Ahj  <Tg-Tw) 

he  =  15.480  P'8  LC^Sp  3 ' 8  B-’2  Tg“'53 

Cj^  -  6.  18  exhaust  stroke 

2.28  combustion-expansion  stroke 

This  study  investigated  heat  transfer  in  engine 
operating  conditions  significantly  below  the  range  for  which 
published  correlations  exist.  Theoretical  work  and  actual 
measurement  of  an  engine  prototype  remain  to  be  completed  to 
prove  the  validity  of  significant  extrapol ati on  of 
reciprocating  engine  heat  transfer  correl ati ons.  This  study 
paints  out  three  specific  areas  requiring  future  study 
concerning  reci procati ng  engine  heat  transfer  prediction: 


lOO 


Quantitative  prediction  of  -flow  pattern  magnitude 
and  distribution  as  a  function  o-f  crank  angle  and 
1 ocati on . 


1 . 


Quantitative  prediction  o-f  soot  concentration, 
apparent  -flame  emissivity,  and  radiant 
temperature  as  a  function  of  operational 
variables  in  low  pressure  diffusion  flame 
combustion . 

Quantitative  effect  of  appreciable  fuel 
injection  and  ignition  delay  in  the  expansion 
stroke  of  low  compression  ratio  reci procati ng 
engines  an  incandescent  particle  concentration 
and  luminous  flame  radiation. 


i. 
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APPENDIX  A 


The  purpose  of  this  appendix  is  to  present  preliminary 
design  parameters  and  the  computer-generated  thermodynami c 
■firing  cycle  -for  the  New  Cycle  Engine  expansion  cylinder  on 
which  the  heat  transfer  calculations  are  based  C1J: 


Preliminary  Design  Parameters;  expansion  cylinder 


1.  Mean  piston  speed,  3p  =  2000  ft/min  <33.33  ft/sec) 

2.  Square  conf iguration;  stroke  =  bore  diameter  (L=B) 

3.  Neglect  clearance  volume;  Yc  =  0 

4.  Neglect  cylinder  sleeve  clearance 

5.  Regenerator  port  exit  cross-secti on  area,  Ap  = 
0.80  ft°  <l/3rd  cylinder  circumference  x  .75  inch 
depth ) 

6.  Compression  cylinder  supercharge  pressure  ratio=l; 
natural  aspirating 

7.  Quiescent  fuel  injection  at  nominal  25°  crank 
angle 

8.  Ignition  start  at  nominal  45°  crank  angle 

9.  Specific  heat  release  per  cycle,  £heat  =  500 

BTU/lbm  charge. 

10.  Specific  fuel  consumption,  SFC  =  .284298  lbm/HP-hr 

11.  Fuel  base  =  nominal  CnH2n  medium  speed  diesel  fuel 

12.  Full  load  equivalence  ratio,  $  =  .75 

13.  Full  load  air/fuel  mass  ratio  =  20 

14.  Mean  heat  transfer  surface  wall  temperature, 

Tw  =  900° R 

15.  Specific  gas  constant  for  air  charge  and 

combustion  products,  R  =  53.90  lbf-f t/lbm-°R 


16.  Air  charge  Prandtl  number,  Pr  =  .70 


TABLE  A-l 


New  Cycle  Engine  expansion  cylinder  computer-generated 
thermodynamic  -firing  cycle  model  as  a  -function  o+  crank 
angle;  5°  increments  Cll 


crank  angle 

Pressure 

mean  7  temp. 

Volume 

0° 

P,  atm 

T  (  0  R 

V,  ttJ 

EXPANSION- INDUCT  ION 

CONDITIONS 

5 

2. 49823 

2149.04 

3 . 406  7 E 

10 

3. 1201 

2159.25 

.  013588: 

15 

3.91693 

2195.73 

. 030424 

20 

4.83479 

nnnc;  n  / 

■  2.0 

.  0537141 

25 

5.67408 

2229. 71 

.  083 1 85: 

30 

6. 10171 

2194.29 

.  118493 

35 

5. 88219 

2116.61 

.  159226 

40 

5. 12078 

2010. 57 

. 204916 

45 

4. 14883 

1896.6 

. 255045 

EXPANSION-COMBUSTION 

CONDITIONS 

50 

4.24684 

2251.6 

. 309053 

55 

4. 09172 

2571.21 

. 366348 

60 

3. 92401 

2864. 13 

.426316 

65 

3. 76064 

3136.02 

. 488331 

70 

3.26107 

3064 . 1 5 

.551767 

75 

2.88976 

2998.73 

. 616005 

80 

2.58662 

2939. 1 1 

■  680446 

85 

2.33694 

2884.75 

. 744518 

90 

2. 12966 

2835. 18 

. 887682 

95 

1 . 95638 

2789.98 

. 869441 

100 

1.81064 

2748. 78 

. 929342 

105 

1.68741 

2711.25 

. 98698 

110 

1.58277 

2677.09 

1 . 042 

115 

1 . 49361 

2646. 03 

1 . 09408 

120 

1.41745 

2617.83 

1.  14298 

125 

1 . 35232 

2592.29 

1.  18847 

130 

1.29661 

2569.22 

1 . 23038 

135 

1 . 24904 

2548.44 

1.26857 

140 

1 . 20857 

aJa.7  ■  Ox. 

1 . 3029 1 

145 

1. 17433 

2513.21 

1 . 33334 

150 

1. 14565 

2498.51 

1 . 3598 

155 

1. 12197 

2485.62 

1 . 38223 

160 

1. 10286 

2474. 46 

1 . 40061 

165 

1 . 08796 

2464.96 

1 .41492 

170 

1.07702 

2457.07 

1.42515 

175 

1 . 06986 

2450. 74 

1.43129 

BDC 

1.06636 

2445.96 

1 . 43333 
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TABLE  A— 1  (CONT'D) 


EXHAUST  CONDITIONS 


P ,  atm 


v,  +  tJ 


1 


BDC 

185 

190 

195 

200 

205 

210 

215 

220 

225 

230 

X- 

240 

245 

250 


1 . 05 
1  -  05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 
1 . 05 


2439. 06 

2434. 63 
2430.7 
2426. 81 
2422. 95 

2419.13 
2415. 36 

2411.64 
2407.99 
2404 . 4 1 
2400. 9 
2397.47 

2394. 14 
2390.89 
2387.75 


1 . 43333 
1 . 43 1 29 
1 . 42515 
1 .41492 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 


o  o  o  o  *—  •-*  to  f.j  (  j  w  u  o-  o  m  oo  co  ~o  -o . . . 

to  tn  co  «-*  tn  o  ui  c  a  to  m  tn  •-  co  o  o-  to  ro  o  o  ►*  ^  to  to  to  to  to  to  -6 
to  O  to  CO  00  03  -6  tn  <3  Ct  O'  00  *-  O'  O  -6  '0  tO  0)  W  -6  -O  -6  03  to  O'  O  tO  tn  CO  O 

u  '  i  —  tn  to  o  o  to  to  to  "j  o  -t>  tn  ^o  £>  to  03*  to  -l  to  ro  o  ro  to  to  o  to  o 

oo  to  to  oo  o  1 1  to  tn  t>  o  to  -6  ro  to  ct  to  tn  tn  03  o  >-*  -u  to  tn  03  to  ro  to  ct 

to  03  to  tn  ro  03  03  03  to  co  tn  to  -6  ct  ro  ro  tn  *-*  03  03  ro  to  tn  to  >-> 


APPENDIX  B 


The  purpose  of  this  Appendix  is  to  present  the 
derivation  of  several  geometric  relations  for  a  simple 
geometry  reel procating  cylinder.  Figure  B-l  illustrates  a 
standard  cylinder  with  defined  piston  linkage  symoology. 

1 .  Instantaneous  piston  speed,  Sp 

by  the  law  of  cosines: 

o  o  r» 

1  ^  =  s“  +  a~  +  2ascos0 

s^  -  (2acosfcOs  +  (a^-l  ‘‘O  =  0 

by  the  quadratic  formula,  solving  for  s: 

-<-2acos©)  ±  C(2acos0)^  -  4  (a^-l  r*")  3  1  /  ^ 
s  =  - 


s  =  acos©  ±  +  a^(cos^-l) 

substituting  in  the  tr i gonometr i c  identity,  cos^y-l=si , 
and  taking  the  positive  real  root: 

s  =  acosO  +  C 1 r2  -  a2sin20 3 1 /2  (B-l) 

differentiating  s  with  respect  to  time,  t: 

ds  dO  t  ds 

de  dt  de 


.  aSinecos© 

Sp  =  9C-asin6  -  -  3 

(lri-aisin‘’0) 

substituting  in  the  crank  ratio  identity,  Rc  =  lr/a,  and 
simpl if ying: 

#  COS0 

Sp  =  -a€>sin6Cl  +  -  3 

r-  —  ^  .  o 

<R=  -si  n^0)  1 


substituting  in  the  &  identity,  6  =  2irNRpM/60: 
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Sp 


-2aNRPM 


60 


(2+0) 


Sp  = 


4aNRPM 


60 


Substituting  in  the  crank  radius  identity,  a  =  L/j 


^lnrpm 

Sp  =  - 

60 

3.  Instantaneous  Cylinder  Volume,  0: 
Neglecting  clearance  volume,  Vc ,  and  sleeve 

V  =  !E&2  Le 

4 


(B— 3) 


clearance: 

( B— 4 ) 


4.  Instantaneous  Heat  Transfer  Area,  ahj! 

Neglecting  clearance  volume,  Vc ,  and  sleeve  clearance: 

aHT  =  2<!l£2)  +  "BL0  (B-5) 

4 

where:  L^=lr+a-s 

s  from  equation  <B-1) 

5.  Calculation  of  New  Cycle  Engine  Expansion  Cylinder 
Preliminary  Design  Parameters: 

B,  L,  1,  a,  R ,  Nrpj^ ,  max  Sp 

From  Appendix  A,  the  following  preliminary  design  conditions 
are  known: 

a.  Clearance  volume,  Vc  =  0  :  L^=lr+a-s 

b.  Sleeve  clearance  =  0  :  B  =  bore  and  piston  diameter 

c.  Square  configuration  :  B  =  L 

d.  Sp  =  33 . 33  ft /sec 

e.  VBDC  =  1.43333  ft3 

f.  v90<*  *  0.80770  ft3 


1 1 1 


At  VBDC’  Lh  =  L  =  B: 

VBDC  =  1  •  43333  -ft3  =  itB2  L0  =  frB3  (B-4) 

4  4 

B  =  L  =  C  4(1.45333)  3 1/3 

IT 

B  =  L  =  1.22204  -ft 

From  the  crank  radius  identity,  a  =  L/2: 
a  =  0. 61102  -ft 

From  equations  (B-l),  <B-4> ,  and  the  geometric  relationship, 

L-0  =  lr  +  a-  sata  given  crank  angle  9: 

4V 

—  s  l  +  a  -  CacosO  +  ( 1  ..^-a^si  n^O)  1  * ^  3 

y  r  r 

ttB 

which  at  known  V,  B,  and  a  values  at  9  =  90°  simplifies  to: 

4 C. 8077) 

- -  =  1 r+. 61 102— C . 61 102cos90°  +  < 1  2-. 61 1022si n290° ) 1 7 

rr(  1.22204)  2 


lr  -  .077615  *  (lr2-. 373346) 1/2 

lr2  -  . 155231 r  +  .006024  =  lr2  -  .373346 

lr  =  2.4440  -ft 

From  the  crank  ratio  identity,  Rc  =  lr/a: 


R_  =  4.0 


From  equation  (B-3) ,  NRpM  =  60  §p/2L: 

60(30.33) 


N 


RPM 


:( 1.22204) 


Nrpm-3..81J§:5 

From  equation  (B-2) ,  recognizing  that  maximum  Sp  occurs  at 
9  =  90° ,  270° : 


Sp 


max 


a2frNRp^si  n90® 
60 


cos90° 


C  1+- 


(Rc  -si n*90° ) I/x 
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FIGURE  3-1 :  Standard  reciprocating  cylinder  and  piston 
linkage  geometry 


APPENDIX  C 


The  purpose  of  this  appendix  is  to  present  formulas  for 
calculation  of  essential  thermal  properties  of  principle 
gases  as  a  function  of  mean  gas  temperature,  Tg,  and  lean 
hydrocarbon-ai  r  combustion  products,  (-'nH2n  fuels,  as  a 
function  of  Tg,  pressure,  P,  and  equivalence  ratio,  §. 

1 .  Principle  gases 

TABLE  C-l 

Coefficients  for  calculation  of  C  and  >*.  of  principle  gases 
as  a  function  of  Tg(°R)  C63 

Valid  for  Tg  >  450° R 


Substance 

AMW 

Y0 

io4y1 

108Yo 

Jim 

107hq 

ai  r 

29 

.224 

.2617 

.  00 

2.  163 

nitrogen  N2 

28 

■  233 

.2478 

.  00 1 2 

2 .101 

carbon  monoxide  CO 

28 

.229 

.  3306 

.2160 

2.053 

carbon  dioxide  CO., 

44 

.  141 

1.3128  -2 

.  4867 

1 . 704 

water  vapor  H-,0 

18 

.479 

-.4483  4 

.  8951 

1 . 499 

AMW  =  Atomic  Mass  Weight 

>*• 

*oT<3 

.  645 

(C-l ) 

lbm/f t-sec 

Cp 

BTU/lbm-®R 

ll 

< 

r. 

+  Y  ^  Tg  + 

Y2(Tg)2 

(C— 2) 

2.  Lean  Hydrocarbon 

-Air 

Combustion  products; 

— £n— 2n- 

fuels 

Valid  for  900® R  i  Tg  i  7200® R  and  f  *=  1.0: 


1 . 470x 10~7 (Tg) ' 7 

)*.  =  -  ( C-3 ) 

lbm/f t-sec  1+.027# 
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*  «  >1 


F'r  =  0.05  +  4 .  2  ( J‘—  1 )  -  6.7<3‘-l)2  (C-4> 

3'(Tg)  s  Figure  C-2 

Valid  for  900° R  ^  Tg  *  3600° R,  P  =  1  atm,  and  #  =  .  75 
(see  Figure  C-2) : 

C  =  .247  +  . 243 x 10~4Tg  (C-5) 

BTU/lbm-0R 
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500  1000  1500  2000  2500 

(900)  (1800)  _  (2700)  „  (3600)  (4500) 

To  :  K  (  R) 


FIGURE  C-1 :  Ratio  of  specific  heats,  X,  of  lean  hydrocarbon- 
air  combustion  products  as  a  function  of  ?,  Tg, 
and  $ ;  fuel  CnH2n 


FIGURE  C-2:  C  of  lean  hydrocarbon-air  combustion  products 
as  a  function  of  P,  7g,  and  $  »  fuel  CnH2n 
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